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O 1 abstract 

O ' 

We study equilibrium solutions for the problem 



in 



c/3 

43 



u t = e 2 u xx — u 3 + \u + cos x 
u x (0,t)=u x (l,t) = 0. 



Using a shooting method we find solutions for all non-zero e. For small e we add to the solutions 
found by previous authors, especially Angennent, Mallet- Paret and Peletier, and by Hale and 
Sakamoto, and also give new elementary ode proofs of their results. Among the new results is 
the existence of internal layer-type solutions. Considering the ode satisfied by equilibria, but 
on an infinite interval, we obtain chaos results for A > Ao = 7^73 and < £ < |. We also consider 
the problem of bifurcation of solutions as A increases from 0. 



^ 2 Introduction 

This is the first of a series of papers studying the existence of bounded solutions of the equation 

O 

e 2 u" = u 3 - Am + cost, (i; 

where A and e are positive parameters. This equation is a standard model in the theory of 
nonlinear oscillations, one of several which have been called a forced Duffing equation in the 
literature. With the signs shown it is often referred to as the equation of a "soft spring." A 
comprehensive reference to the early theory of (U) is [ |NM|| , which gives a detailed account of the 



results obtained by classical perturbation methods, such as averaging or multi-scale techniques. 
More recent efforts have used dynamical systems concepts to establish results about chaotic 
behavior of one sort or another. An important reference is by Angenent, Mallet-Paret, and 
Peletier |[AMPP|1, who studied stable steady states for a react ion- diffusion equation 



u T = e 2 u xx + f (x,u) (2) 

with boundary conditions 

u x (0,r) = u x (L,t) = (3) 

for a class of functions / which were cubic in the state variable u. Equation (|I|) is obtained 
from (0) (with / (x, u) = u 3 — \u + cos x) by setting x = t and assuming that u is independent 
of r. While the specific form (HD is not mentioned in |[AMPP 1, the methods there are easily 
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seen to apply and to prove important results about the existence of stable steady-states of the 
problem (0) — (|3|) with this function /. Compactness requirements, however, seem to make it 
more difficult to study the problem on an infinite interval, and thereby get "chaos" , using infinite 
dimensional methods. 

In this paper we introduce a shooting technique, which we have not seen in this form elsewhere, 
which is the basis of our approach. We use this technique to obtain new periodic solutions for ([J) 
over a range of e which is not "small" . We show that one can give rigorous results about a weak 
form of "chaos" over this larger range of e. We also introduce the problem of determining how 
these solutions arise as A varies, since for A < 0, (1.1) will be shown to have a unique bounded 
solution. Then, letting e return to its traditional role as a small parameter, we reproduce and 
extend the results in [ |AMPP|| as they apply to this equation, using elementary ode methods. 



We study the bifurcation problem in A in more detail, and we begin our study of a further class 
of bounded solutions which we can find for small e. 

However, a detailed study of this new class of solutions has proven rather lengthy, so we will 
leave its elaboration to a later paper. The goal of that paper will be to obtain a general chaos 
result which includes all of the solutions discussed in each paper, in some sense. It is anticipated 
that there will be still a third paper devoted to this topic, in which we use techniques from 
HM2|| to study the stability of the new solutions, as equilibria of (0) — (H) . This will lead to a 



consideration of the Morse index of these solutions, and their role as part of the global attractor 
for (0) — (|3|) . Future work will include extension of the results to more general equations. 

In most of this paper we will consider the ode (H) for its own sake. The literature is large, 
and it is surprising that, as far as we have been able to determine, many of the solutions we 
will find have not been described before. Some of these solutions are periodic and some are 
"chaotic" in a sense to be defined below. Those which are periodic are unstable as steady-states 
of (fj) — (Q) and therefore apparently not found by the pde techniques of ||AMPP|| . However recent 



work of K. Nakashima on a similar problem uses continuation methods which may yield results 
similar to some of ours for this problem, such as the existence of these new periodic solutions 
JNAK]. There have been many other papers on the problem (|2[) — (^) with cubic nonlinearities, 



and some of these also consider unstable solutions. Two recent ones, which cite other related 
work, are by Hale and Salazar ||HS1|1 , |[HS2|1 . An earlier paper by Kurland gives many oscillatory 



solutions which are similar to some of those we find in section 4.2, but for a different class of 



nonlinear functions f(x,u) [|KUR|| . The use of formal asymptotic analysis on similar problems 
has been studied by Ockendon, Ockendon and Johnson [ UUJ |, Norbury and Yeh | N Y| | and Mays 
and Norbury |Mg . 



All of these papers use methods which appear to quite different from ours, and study, for the 
most part, equations different from ([!]). We have not investigated whether these methods, many 
using infinite dimensional functional analysis or sophisticated topology, apply to ([]]). Most do 
not mention equations of the general form of ours, (as in (|5|) below), ||AMPP|| being a notable 



exception. Most do not study problems on an infinite interval, so that chaos is not considered. 
We have seen no other work on the bifurcation problem considered in section 4.5. 

The proof of a weak form of chaos which we give seems to us to be very simple. No analysis 
is required beyond a simple phase plane argument and the continuity of solutions with respect 
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to initial conditions. See Theorem [12] and also the first part of section [4.4| . This is for small 



e. For larger e a few estimates are needed to verify the hypotheses of Theorem 12. These also 



appear in section \LA 



To relate our results to those of traditional nonlinear oscillation theory, we recall that in HNMfl 
the undamped form of the equation is written as 

u + ujqU = e (u 3 + k cos ((w + ea) t)) . (4) 

This form is chosen to study "near-resonance" phenomena, which were the main interest of much 
previous work. Resonance, or near-resonance, will play no role in our approach, since we are 
not studying a small perturbation of a linear problem. A rescaling of (f|) to put the equation in 
the form (|I|) shows that e small in (^) corresponds to A large in (|T|) , though the analogy is not 
exact because this rescaling also results in a small amplitude high frequency forcing. 

The paper is organized partly according to the restrictions place on e. Propositions III, 0, 
and Theorem |] in section [| are for any e > 0. The results in Section 3 are for a range of e 
which can be stated explicitly, while the results in section 4 are for "sufficiently small" e. A 
more detailed outline of the paper can be found in section 5. 



3 Results for all positive s 

At this initial stage it is easy, and we believe of some interest, to study a more general class of 
equations. So to start with, consider the equation 

e 2 u" = u 3 - Xu + g (t) (5) 

where g is any continuous bounded function on [0, oo), with initial conditions 



u(0) = a, u'{0) = 0. (6) 

We denote the unique solution by u a . Our initial interest for a general g is to find solutions which 
are bounded on [0, oo). Fixing e > 0, this can be viewed as a bifurcation problem in A. There 
is a "main branch" of bounded solutions which exist for all A, as given in the following result. 



Proposition 1 For any e > and any A, (§D has a solution satisfying (||) which is bounded 
on [0, oo). 

Proof. Choose b > so large that if \u\ > b then u 3 — Am + g (t) ^ 0, for any t. Then, 
u" > when u > b and u" < when u < —b. We define two subsets of the a axis: 

A = {a | u a (x) > b for some x > 0} 
B = {a | u a (x) < —b for some x > 0} . 
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These sets are clearly non-empty (e.g. b+ 1 G A) and open. They are disjoint because u a cannot 
have a local maximum in the region u > b or a local minimum in the region u < —b. Hence there 
is an a* ^ A U B, and this corresponds to a solution which is bounded on [0, oo). This proves 
Proposition [I]. ■ 

It should be remarked that for any even function g (•) , such as cosine, the solutions u a are 
even and so u a * is bounded on the whole real line. 

Proposition 2 For any A < there is only one solution of (§) — (§) which is bounded on 
[0,oo). 

Proof. Suppose that there are two bounded solutions, say U\ and and let v — u% — u 2 - 
Then 



e 2 ~" 



v" = (u\ + U\U% + w 2 , — A) V, V (0) = 0, 

and since u\ + U\U 2 + u\ is positive definite and A < 0, v cannot be bounded on [0, oo). This 
proves the result. ■ 

Remark 3 A similar argument shows that for a general continuous bounded function g, when 
A < there is exactly one solution of (0) which is bounded on (— oo, oo) . However the existence 
part of this result is a little longer, requiring a two-parameter "shooting" argument and some 
more complicated topology. 

We now consider what happens as A increases from zero. Consider the function p\ (u) = 
u 3 — Am. As A crosses zero, p develops, by a "pitchfork bifurcation", two new roots (besides 

u — 0), a local maximum at u — ~\pi an d a local minimum at u — y|. As A increases further, 
the value of p\ at these maxima and minima will eventually exceed the maximum of \g (t) \ . Let 

Aq = sup { \ \ \g (t)\ > p\ I — \ ) for some t 



Proposition 4 If X > Ao, then for any e > there is a unique bounded solution with u > 0, 
3m 2 — A > , and another bounded solution, also unique, with u < 0, 3m 2 — A > 0. 



Proof. Existence of these solutions is a simple application of the shooting technique from 
Proposition |]. It is easy to show that if A > Ao. then at u = where u 3 — Xu has a local 

minimum, u" < 0. We consider solutions u a with \ / 1 < a < b. For a in this region and close to 



|, u a decreases and crosses y | transversely. For a in this region and close to b, u a increases 
and crosses u = b, again transversely. It easily follows that there is a solution u a such that 
\/| < u a (t) < b for all t > 0. Similarly, a second solution is found with —b<u a < — yf • 
Uniqueness of both of these solutions follows as in Proposition |2|. ■ 
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Remark 5 We believe that there should be a third bounded solution as well, but we don't have a 
proof for a general continuous bounded function g. There appears to be a significant difference 
between boundary value problems on finite intervals, such as those in j\AMPP\l , and the semi- 
infinite interval problem discussed here. It is easy to show, for a bounded g and any given 
L > 0, that for large enough A there are at least three solutions of (^) — (§) such that u' (L) = 0. 
But it does not appear to us to be so easy to find three bounded solutions on [0, oo). "Generically" 
there should be three, perhaps by some sort of degree theory argument, but the non- compactness 
of [0, oo) seems to make the use of such arguments more challenging. 



So now we specialize to the equation of principle interest, namely ([[]) . It is useful to consider 
the curve in the (t, u) plane defined implicitly by the equation u" = 0; that is , by / (t, u) = 
u 3 — \u + cost = 0. This curve has one component for A < An = \ and three components for 
A > An. Here arethe graphs for < A < A 0j A = A , and A > A : 




Figure 1. 

In all three cases, the curve crosses u = at odd multiples of |. For A > An, let U (t) > 
Uq (t) >U_(t) be the three solutions of u 3 — \u + cos (t) = 0. 



Theorem 6 When g (t) = cost, letu\ andu$ be the solutions found in Proposition^, withu\ < 
and m 5 > 0. Then these solutions are 2tt -periodic, u[ > in (0,7r) and U_(0) < U\ < U_(tt) in 
[0, 7r] . Also, W5 (t) = —Mi (t + 7r) and u\ and u$ are the minimal and maximal bounded solutions 
of ([I]) on (— oo, oo) , in the sense that every other bounded solution lies between these two. Also 
there is a third 2tt periodic solution, u p (= u a A , which satisfies u (|) = and u' > on (0, |). 
A unique solution with these properties exists for all A (and so for A < this is the solution 
found in Proposition^). 

Proof. In proving that the solutions of Proposition f| are periodic, we give a second proof of 
their existence, for the case g (t) = cost. Let a = U(0). First observe that u'^ Q (0) = u^ (0) = 

while e 2 Uao(0) = —1. Therefore u ao (t) decreases monotonically to — oo. If a < a < U (0), then 
w"(0) > and so u a initially increases. However, for sufficiently small positive values of a — an, 
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u' a has a first zero, which we denote by t\(a). This means, in turn, that u" a must have a first 
zero at t\ (a) G (0, t\ (a)). For small a — a , t\ (a) < |. Further, e 2 u'^ {t\ (a)) < 0. From 

£ 2 u"" (t) = (3m 2 - A) u" + 6mm' 2 - cost 

we conclude that for small a — a > 0, v!' a < on (ti (a) , ti (a)}. Hence, u (t 1 (a)) < U_ (t 1 (a)) . 

Clearly, ti(a) approaches as a tends to a from above. Extend the function tx(-) to larger a 
by continuity, as long as possible. That is, ti(-) is the continuous function such that u' Q (ti(a)) = 
and, for a sufficiently close to «o, t\{a) is the first positive zero of u' a . Then t\ (a) remains the 
first zero of u a until either (i) u'^(ti(a)) = (since we can then no longer use the implicit 
function theorem to solve u' a (t) = for t) or (ii) there exists a first to £ (0,£i (a)) such that 
u' a (to) = and u'^(to) = (since then t\ (a) ceases to be the first positive zero of u' a ). But 
if (i) occurs at some first a\ > a , and t\ («i) is the first zero of u' with ti («i) < tc, then 
M ai (*i ( a i)) = ~ sinti (ai) < 0. This implies that u' ai = to the left of ti («i), a contradiction. 
If (ii) occurs we get the same contradiction at t (a) . If (i) occurs at a\ and t\ («i) = ir, then 
M ai (*i ( a i)) = an d £ 2 «a'i' (*i ( a i)) = 1- This means that v!' ax > on either side of t\ [a\) = ir, 
and t\ («i) could not be the first zero of u'. Hence ti (a) is continuous and remains the first 
zero of u' a as long as t\ (a) < n. 

On the other hand, if a = U_{^)i then the first maximum of u a is larger than Uo(7r). Therefore, 
ti(a) must vary continuously as a increases from cto until it eventually takes on the value 7r for 
some a G {<yo,U_(n))- For any such a, u a gives a periodic solution of (P with the properties that 
u' a (t) > for t G (0,7r) and C/(0) < w Q (t) < C/(7r) for t G [0,7r]. From Proposition |, it follows 
that such an a is unique and the corresponding solution u a coincides with the unique bounded 
solution with u < and 3m 2 — A < in Proposition f|. The reflection — u a (n + 1) of this periodic 
solution gives a second periodic solution, which corresponds to the unique bounded solution with 
u > and 3m 2 — A > in Proposition |j. 

The proof that Ui and M5 are the minimal and maximal bounded solutions of ([ft) on (—00, 00) 
follows from that of Proposition 

To obtain the third periodic solution, u p , again let «o — LL (0) , the smallest root of u 3 — Xu + 
1 = 0. We saw that if a < a then u 3 — Xu + cost < for all u < a, and all t, and u a decreases 
monotonically to —00, (in finite time). If a = a , so that u" (0) = 0, then u'" (0) = 0, u"" (0) < 
and again u decreases below a and tends to —00. On the other hand, if a < is sufficiently 
small in magnitude then e 2 u" (0) > | and it is easy to show that u increases to cross zero at some 
first to — ^0 ( a ) > with u' > on (0, to («)]• Lowering a from 0, the function t (a) is continuous 
as long as u' a (to {a)) > 0. If u' a (t (a)) = then u" (t (a)) < since u cannot have a local 
minimum at to (a) . This implies that to (a) > |. Since to (a) is not defined for a < «o, it must 
increase beyond | as a decreases from zero. Thus there is some first (largest) a = a p < with 
( a p) — f • Let u p = u ap . An easy symmetry argument shows that u p (| + t) = — u p (| — t) for 
all t, showing that u p is 2ir— periodic, with a maximum at 7r. 

Suppose, therefore, that u = u p satisfies u' (r) = for some first r G (0, n) . By the anti- 
symmetry of u p around | we can assume that r < |. Then (r) < 0. If u" (r) = then 
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e 2 u"' (t) = — sin r < so u" becomes negative. In either case u' < on some interval to the 
right of r. 

From the graphs of u" = it is clear that if u (0) < 0, u"(0) > 0, u < on (0, f ) and vl = 
before t = | , then after that point, it" < at least until t = ir. This implies that u' p > on 

(0,7T). 

The uniqueness of u p follows from Lemma [7] below. 

Lemma 7 Assume that a% < (%2 < 0. If < T < n/2 and u a2 < in [0,T], i/ien w ai < u a2 in 
[0,T]. 

Proof. First, using the Sturm-Liouville technique we get 

e 2 (u' ai u a2 - u ai u' a2 )(t) = / u ai u a2 (u 2 ai - u 2 a2 ) ds + / (u a2 - u ai ) cost ds. (7) 

Assume that the lemma is false. Then, there is a first t = i G (0, 7r/2] such that M a2 < in 
[0,t] and u ai = u a2 at t. Then, u'(t) > u'(t). Evaluating (0) at t — t shows that the left side 
is nonpositive and the right side is positive. This contradiction proves the lemma and completes 
the proof of Theorem []. ■ ■ 

Below we will want a slight variant of Lemma 0, with a similar proof which we omit. 

Lemma 8 Let u and U be two solutions of with u'(tt) = U'(n) = 0. If n > T > n/2 and 
u > in [T, 7r] , then u < U in [T, n] . 

Corollary 9 Assume that u is a periodic solution of (fjj with u'(0) = u'iji) = 0. Ifu(0) < u p (0), 
then u{t) < Up{t) for all t G (— oo, oo). 

Proof. From Lemma [7] it suffices to show that u(t) < u p (t) for t G (tc/2,tt]. If u(n) > u p (tt), 
then Lemma ^| implies u(t) > u p (t) for all t G [| , 7r] so that u has a jump at |, a contradiction. 
Hence we have u(7r) < u p (tt). Assume now that the corollary is false. Since u' p > in (0, 7r), 
there is a t G (|, tt) such that = and u(t) > u p {t). The Sturm-Liouville technique shows 
that 

e 2 [u \t)u p {t) — u(t)u' p (t)] = —u(t)u' p (t) + J [uu p {u 2 — u 2 ) + (u p — u) coss] ds. 

If u(t) = u p (t) for some largest t G [f,t] then < in (t, t) a contradiction. We again get 
that w has a jump at t — ~ and this proves the corollary. ■ 

We have now shown that as A increases from 0, new periodic solutions appear. We can 
prove one thing about the initial bifurcation of new solutions for any positive e. Let 

\b = sup {A | There is only one solution, u p , with u' (0) = 0, u (tt) = 0} . 
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Theorem 10 For any e > 0, A b < A 



of solutions of u" = 0, shown in Figure 1-b. Starting from a = — J \ we consider the solution 



Proof. Suppose that e > and let A = A . It will be helpful to consider the graph of the set 

3 

u a as a is lowered. Let t 2 («) denote the first intersection of this solution with u = —\J\- 

Then near a = — <y/|, t 2 (a) is defined and continuous. Further, u' a > on (0,t 2 («)] as long 

as t 2 (oc) < 7T. The solution cannot be tangent to the line w = — y| until we reach a value of a, 
say ai, where £2 («) = 7r. Suppose that u' a (n) = 0. Since A = Ao it follows that u'^ (it) = 

and we easily calculate that u'^n) = while w™ (tt) > 0. But this implies that u ai = ~\pi 
at some earlier point, a contradiction. Therefore u' a (re) > 0. A similar argument shows that 

there is an a 2 > ■^/| such that u a , 2 {n) = and u' a2 {ji) < 0. If a > b then u' a (ir) > (or u a (n) 

is not defined, which is easily handled), while if a < — b then u' a (ir) < 0. Therefore a simple 
continuity argument shows that there are three values of a such that u'(tt) = 0, and this give 
three 27r-periodic solutions of [l]. 

Further, this argument extends by continuity to values of A close to Ao- That is, keeping ot\ 
and «2 fixed, we still have u' (n) > 0, u' a2 (n) < 0, u' b (ir) > or undefined, and u'_ b (n) < or 
undefined. This proves Theorem ■ 

4 Results for "large" e 
4.1 A simple proof of chaos 

In this section we show how the shooting method can give a global result about chaos for 
this equation. Our aim is for a concise statement implying the existence of uncountably many 
bounded solutions and a natural map between a family of symbol sequences and a set of bounded 
solutions, where moreover, we can get specific estimates on ranges of A and e which support this 
chaotic behavior. This will mean that we do not prove uniqueness of the solution corresponding 
to a particular sequence, and hence we can only prove a weaker sensitivity to initial conditions 
then one obtains from some asymptotic methods. In a later section, in which e is taken as 
sufficiently small, more precise results will be obtained. 

Our results in the next few sections have one main hypothesis, which we state now. 



Condition 11 There is an a G (cto, — suc ^ that u & (t) increases monotonically as t 

increases from 0, and crosses u = b before t — | . 

Recall that b is positive and satisfies b 3 — Xb — 1 > 0. An easy phase plane argument (given 



in section |4.4j) shows that if A > Aq, then Condition ITT] is satisfied for sufficiently small e. Later 



we will show, with a bit more work, that Condition 11 is satisfied in a specific range of e and A. 



8 



We use a as given in Condition |TT] to define a family {lik} of "special" solutions as follows: 

w k {t) = {-l) k u a {t-k7i). (8) 

Thus, each Wk is a translation, and for odd k a reflection, of the solution u s . It is easily seen 
that (— l) k Wk has its global minimum at kn and on each side of this point increases monotonically 
to u = b before \t — kn\ = |. 

Let Wk denote the restriction of Wj~ to the interval [sf., Sk] in which \wk\ < b. In the following 
figure we plot several of the functions Wk- 



Figure 2 
Graphs of w ,wx,W2,w 3 , and u> 4 



We can now state our result. 



Theorem 12 Assume that A > Ao and that Condition [TJ holds. Let X = {o~fc} fc=1 x be any 
sequence of positive integers such that — o~ k > 2. Then there is an «s G (ao, a) such that 
u as exists on (— oo, oo) , and on (0, oo) the graph of u as intersects the graph of each w Uk , but 



does not intersect the graph of any Wj with j E. 



Remark 13 We allow S to be finite, or even empty. 

Proof. We define two functions, /+ and /_, as follows: 

„ / ,\ _ J Wk (t) if k is even and Sk < t < Sk 
+ \ b otherwise 

„ /,\ _ J Wk (t) if k is odd and Sk < t < Sk 
\ —b otherwise 

A key observation is that no solution can intersect the curve /_ tangentially from above, or the 
curve /+ tangentially from below. This is because u" < when u = —b, u" > when u = b, and 
no two distinct solution graphs can be tangent to each other. Also, no solution can be tangent 
to either of the lines u = ±b from between these two lines. 
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Further, we will need two other functions, g + and g_, where g~ < < g + . These are defined 

by 

9+ (t) = max | /_ (t) | , g- (t) = min | f+ (*) J • 

As in the figure below, the function / + can be described as the function whose graph is the line 
u = b except at downward bumps when this line meets a w k with k even. Also, the graph of 

g + is the line u = y^| except for upward bumps when this line meets the graph of w k for some 
odd k. The functions /_ and g_ are reflections and translations of / + and g + . Recall that in 
Condition 111] we required that |a| > \ \. No solution can be tangent to g + from above or to 



g>_ from below. Also, no solution can be tangent to u — y | from above or to u = — y ~ from 
below. 




Figure 3 



Graphs of /+ and g + 

Assuming that the sequence £ is chosen, we shall define a sequence of closed intervals {I k } 
with the following properties: (See Figure 4). 

(z) For each k, a £ I k implies that u a intersects w ai , ...,w CTk and no other Wj with j < a k . 
(ii) Let Ik = [Pk, Oik] , an d for a G I k let t k (a) denote the first t such that u a (t) = w Uk (t) . 
Then either t k (a k ) = s ak and t k ((3 k ) = S ak , or t k (f3 k ) = s ak and t k (a k ) = S ak . 
(Hi) I k C 4_i. 

Since the intervals will be closed, their intersection will be non-empty, and a point as G nj^jlfc 
will have the properties stated in the Theorem. We do not prove that the intersection is only 
one point. We will be able to do this for sufficiently small e. 

We will choose I\ as a subinterval of (—6, a) . Suppose, first, that o\ is even. For a close to 
a, u a intersects f + before t = ~, and therefore before u a could intersect W\. Let t\ (a) denote 
the first intersection of u a with / + . Then t\ (-) is continuous in some maximal subinterval (/3, a] 
of (—6, a] . However, j3 > —b, for solutions u a with a close to —b do not intersect /+ at all. 
(Instead they decrease monotonically to below u = —b and then to — oo.) The nontangency of 
u a with f + implies that 

lim sup t\ ((3) = oo. 
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We do not know if ti (a) is monotone increasing. However, there must be an interval l\ = \J3i, oti] 
such that ti (ai) = s ai ,ti (Pi) = S ai , and s ai < t 1 (a) < S ai in (f3i,ai) , ti (a) < S ai on (Pi, a] . 
As long as ti (a) is continuous, it must remain the first intersection of u a with / + , because of the 
nontangency of u a with / + . Further, while ti (a) is continuous the solution u a cannot intersect 
/_ in (0,ti (a)) , because of non-tangency with /_. A similar argument is used to construct Ii 
if cr 1 is odd. In this case we have = s ai and ti(ai) = S ai . 

We now assume that a (decreasing nested) sequence of closed intervals Ii,...,I n has been 
constructed with properties (i) — (iii) for k = 1 • • -n. We wish to construct I n+ i C /„. We 
will consider the case where a n is even and u n+ i is odd. We will also assume that t n (a n ) = s an 
while t n (j3 n ) = S an . We will show that we can construct [ct n+ i, j3 n +i] so that t n+ i (a n+ i) = S Cn+1 
and t n+ i (Pn+i) = s 

tn+l (Pn+1 



s, 



, or we could make a different choice which would result in t n+ i (a n+ i) = 
Un+1 . (Hence, for any given sequence £ there will be many solutions with 
the properties (i) — (iii) . This does not by itself imply a reduced sensitivity to initial conditions, 
because these solutions may be separated from each other. This is discussed in a separate section 
below. ) 

The assumption that a n is even means that the function w an is a downward pointing spike 
and w an (s an ) = w Cn (S an ) = b, while taking a n+ i odd means that w an+1 is an upward pointing 
spike. 





"A '- 













Figure 4 

Several solutions used in the shooting process, with the scale t — > |, e — .5 
Solutions shown are 11^,11^, w m and uo l , together with wo, wi, W2, and W3. 

To define I n +i, in the case where a n is even, let [i n = sup ja < a n \ u a (t n (a)) = • Then 
t n (a) G (s an ,a n ir) , u a (t n (a)) G (-y^f, for < « < ot n , and u 1 ^ (t n (//„)) < 0. The last 
inequality is true because u^ n does not intersect t(V n _i, an upward pointing spike, and \a\ > 
(If we wished to have t n+ i (a n+i ) = s Un+1 we 



would set /j, n — inf |/3 > (3 n \ u a (t n (a)) = -y/f | • 

In this case t n (a) e (a n n, S n ) , u a (t n (a)) G & ) for (3 n < ot < [i n , and u'^ (t n , fi n ) > 0. 

Slight changes are necessary if a n is odd.) 



Consider values of a near to pi n . Since u Un (t n (/i„)) = y | , s n < t n (jj, n ) < cr n TC, and 
(t n (/i„)) < 0, we can define r n (a) continuously by the equations u a (r n (a)) = g + (r n (a)) , 
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Tn (yUn) = t n (fi n ) . Then r n (■) will be continuous in some maximal interval around p n - However, 
as a increases from p n , r n (a) must tend to infinity, since it is not defined at a = a n , where the 
solution u reaches b at s n = t n (a n ) and never decreases below b after that. 

Let 

9 n = sup {a > ji n \ r n is continuous on \p, n , a] and r n (a) = (a n + 1) 7r} . 

The interval I n +i will be a subinterval of [9 n , a n ] . 

Observe that ug n crosses w an+ i, from above, at its maximum point, at t = (cr n + 1) tc. This 
is also an intersection of ug n with /_. Let the intersection of this solution with /_ be denoted 
by p = p n {On), and extend p n as a function of a continuously for a > 9 n , p n being defined as 
the intersection of u a with /_ , as long as it remains continuous. The solution u a may possibly 
intersect /_ at earlier points, but p n is defined uniquely by requiring that p n (9 n ) = (a n + 1) it. 

We see that as a increases from 9 n the function p n must eventually be undefined, since it is 
not defined at a n ■ So, it must increase (not necessarily monotonically) , and there must be some 
closed subinterval of (9 n ,a n ) in which p n (a) lies in the interval [s cr „ +1 , S , (Tn+1 ] and moves from 
the left end of this interval to the right end (not necessarily monotonically) as a increases. This 
subinterval of (9 n , a n ) is chosen for the interval I n +i = [/3 n +i, &n+i]- We define it unambiguously 
by choosing the subinterval with the given properties which lies nearest to 9 n - The possibility of 
an infinite set of subintervals with these properties is precluded by bounds on the variables for 
a given e. The construction of I n +i shows that it satisfies conditions (i)-(iii). 

A similar construction will give I n +i in the case where a n and a n +i are both even. In this case 
we can obtain t n+ \ (a n +i) = s a „ +1 by decreasing a from a n and observing that the intersection 
of u a with u = f + must tend to oo. before we reach 9 n - The case where a n is odd is also handled 
similarly. At each step we can obtain t n (a n ) = s an if a n is even and t n (a n ) = S an if a n is odd. 



This completes the induction step and the proof of Theorem 12 



4.2 "Kneading" theory 

A one-parameter shooting process as in the proof of Theorem [12] inherently gives an ordering 
of the initial conditions corresponding to different sequences. This sort of ordering is related 
to "kneading theory" |GH|| ■ In earlier work on similar problems formal results about this were 
stated, for symbol sequences on two symbols [|HM1|| , [HM2] , [pTT| , ||HT2|| . Here we will describe 



the theory only for the solutions found in Theorem [12]. Further solutions found below would 



make the description more complicated and so this part of the theory will not be discussed in 
those cases. 



Notice that in the induction procedure used in the proof of Theorem |T^, repeated use was 
made of lowering or raising a and letting the intersections with the sets /+ or /_ tend to infinity. 
There was never an assumption that these intersections moved monotonically. However, the 
intersection will pass through some particular Wk before it first intersects Wk+2, the next spike 
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pointing in the same direction as Wk- The point of intersection might retreat along the t axis 
and intersect Wk again. This would increase the number of solutions corresponding to some 
particular sequence, but does not prevent our giving an order to those chosen as in the proof of 



theorem [12]. Our repeated choice to have t n (a n ) = s an , rather than S an , makes the relation 
between solutions and sequences easier to describe. (If o\ = 1 then we could not make this 
choice at step 1, but could thereafter.) The ordering of solutions corresponding to sequences by 
the algorithm used in the proof is well defined because of the definitions of fi n , 9 n , etc. This 
is despite the fact that we don't know that the nested intervals converge to single points. Any 
choice of points within the limit intervals will give the same ordering. We have the following 
corollary of the proof of Theorem O: 



Theorem 14 Suppose that Si and £ 2 are two sequences of positive integers as in the statement 
of Theorem [7^. Suppose that for some k > 0, 0i(£i) = crj(S 2 ) if 1 < i < k. If ak+ii^i) 
is even and <7fc+i(£ 2 ) is odd, then a(£i) > a(£ 2 ). If <7fc + i(£i) > (Tfc + i(£ 2 ) and both of these 
integers are even, then a(£i) < «(S 2 ). If Ck+ii^i) > crfc+i(S 2 ) and both integers are odd, then 
a(Ei)>a(S 2 ). 



4.3 Further periodic solutions, and symbolic dynamics with five sym- 
bols. 

Earlier we showed that for A > Ao, there are at least three solutions with period 2tt. Numerically, 
however, we find five periodic solutions u±, U5 , with u\ < u 2 < u 3 < m 4 < u 5 Here is a graph 
of these solutions when e = 1, A = 2: 




Figure 5 



Periodic solutions, Ux, w 2 , M3 (= u p ) , w^and M5. Here e = 1, A = 2. 



The solution w 3 is the solution u p which is anti-symmetric around |, found in Theorem || 
The solutions U\ and M5 are the ones found in Proposition |j. These are the same ones found in 
AMPPf , at least for small e, and correspond to stable steady states for (|2|) — (|3|) . The solutions 



ii 2 and 1*4 seem to be new (although their existence is strongly suggested by the existence of three 
stable steady-states). Other solutions are found which "follow" one or another of these solutions 
over long intervals, switching in almost arbitrary ways from one to another. To state a precise 
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theorem, we have to capture what it means to "follow" a solution over an interval in a way 
that sharply distinguishes different patterns. We need to identify disjoint open sets of solutions 
which, over a finite interval, follow the solutions in different orders. We are able to do this 
because of the various non-tangency conditions stated above. 

The following result is stated without reference to ui, or even knowledge that u% and u± 
exist. It will be seen later that existence of u<i and M4 is an easy consequence of the technique 
used to prove the result. 



Theorem 15 Suppose that the hypotheses of Theorem are satisfied. (That is, A > Ao and 
Conditioning holds.) Let {uik} be the family of "spike" solutions defined by equation ([8]). Let 
Q = {uJk} be a sequence chosen from the integers {1,2,3,4,5}. Assume that no 1 or 2 is 
followed immediately by a 4 or 5 ; and no A or 5 is followed immediately by a 1 or 2. (There 
must be an intervening 3.) Then there is a solution u = uq with the following properties: 

(z) If uJk — 1 or 2 then u intersects ui2k-\ o,t two points and does not 
intersect W2k- 

(ii) If u>k = 1 then u < g_ in {{2k — 2) it, 2kix) while if uj^ = 2, then u > g_ somewhere 
in this interval. 

{Hi) If uJk = 3 then u does not intersect either W2k-i or W2k- 
{iv) If uik = 4 or 5 , then u does not intersect W2k-i but does intersect W2k, at two points, 
{v) If UJk = 4 then u falls below g + in {{2k — 1) tt, {2k + 1) tt) , while if uik = 5, then 
u > g + in this interval. 



Proof. The proof is a refinement of the proof of Theorem |12[ The notation is different 



however, for in the statement of Theorem [15] all of the Wk are included, not just those picked 
out by some index set {o-fc} . As before, Sk < will be the endpoints of the graphs of Wk- 
The functions /_ and /+, g_ and g + , are the same as before. Also, t k (a) will denote the first 
intersection of u a with 

Rather than carry out a formal induction, we will show how to construct a solution corre- 
sponding to a sequence beginning with uji = 4 or 5, uj 2 = 3,uj 3 = 1 or 2. Hence we want a solution 
which misses wi, intersects W2, misses W3 and W4, intersects w$, and does not intersect wq. It 
will then be apparent how to deal with the general case. In the notation of Theorem [T2| we are 
starting with o\ = 2, 02 = 5. 



If uji — 5, then we proceed in the same way as in the proof of Theorem [L2[ The induction 
step there (going from n — 1 to n — 2) produced an interval I2 in which the solutions u a do 
not cross g + in [tt, 3tt] . Further, the solutions in this interval will not intersect W3 or W4 and 
will intersect the two endpoints of when a is at the endpoints of 12- The refinement of this 
interval according to whether uj% = 1 or 2 will be similar to the way we handle the case Wi = 4, 
so we turn to that case now. 
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When u\ = 4, we follow the proof of Theorem [L2| up to the definition of fix- However now we 

set 




fix = inf <^a > Pi\ t\ (a) < 2tt and u a (tx (a)) = 

Then for a G [fix, fix] > u a intersects W2 to the right of the first point where u> 2 = \/|. We now let 
Tx (a) be the intersection of u a with g + defined by continuity from the point (tx (fix)) ■ Set 

6x = inf {a > fii\ Tx (a) = 3tt} , 

and 

4>x = sup {a < fix\ tx (a) = 2tt} . 

These are defined because u ai does not intersect g + after tt, and because tx (a) varies between 
Tx (fix) and Sx as a moves from fix to f3\. Then <px G (Pi,Qi) and for a G [</>i,#i) u a intersects 
w 2 and also intersects g + in the interval (it, 3tt) . 

The solution Ufa intersects W2 at its minimum point, at t = 2tc, and from there increases to 
cross b before 4p (This is shown by a comparison with w 2 , as in the proof of Lemma [1].) The 
solution uq 1 intersects w 2 twice, intersects w 3 at its maximum point and then decreases to cross 
—b before -y. Following the same procedure as in the proof of Theorem we lower a from 
Ox and find an interval I2 = [$1,012] C (4>x,9i) such that u a2 (S3) = —b, up, 2 (S3) = —b, and 
for a G (l3 2 , 0L2) , u a intersects g + in (7r, 3tt) and also intersects W3. Other cases, and succeeding 
steps, being similar, this completes the proof of Theorem O. ■ 



We can also use this process to prove the existence of U2 and u^. 



Theorem 16 If\ > \q and Condition\T^ holds, then there are at least five 2n- periodic solutions 
of H], Ux, u 2 , u 3 (= u p ) ,m 4 and u 5 , all with u' (tt) = 0. Further, these solutions are ordered by 

Ux < U2 < < U4 < M5. 

Proof. We already know that ux, U3, and W5, which exist by Proposition f|, are increasing in 
(0, 7r) . We give two proofs of the existence of U2 and U4. The first is a little shorter, while the 
second shows in addition that w 2 and U4 are increasing on (0, tt) . 

Proof 1: We will find a solution u = u a with the property that a < a p , u' (tt) = 0, u (tt) > g-, 
and u intersects Wx- (None of the three solutions u 1: u 3 , and u 5 have all of these properties.) 
We assume in Theorem O that lux — 2. Then in the construction used in the proof of Theorem 



T5| we find fix such that u' > on (0, tx (fix)] and w w (tx (fix)) = 9~- We also find <px such 
that tx (4>x) = tt, and necessarily, u'^ (tx (4>x)) < 0. If we decrease a from <px we come to a first 
(largest) a, say a, where u& first intersects g>_ at t = tt. At that point, u' & (tt) > 0. Somewhere 
between a and 4>i there must be an a with (7r) = 0, u a (tt) > g_. This solution must also 
intersect u>i, and this gives the desired fourth periodic solution (beyond those constructed in 
Proposition |j.) The fifth is obtained by translation and reflection. From Corollary || it follows 
that ux < U2 < W3 < U4 < M5. 
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Here is a second proof: If a = Wi(0), then v(t) = ^ |(t) > and v'{t) > for t G (0, 7r], and 

itj > in (0, 7r). Also, Ui intersects W\ in (0, it), and it follows that for a > Ui (0) and sufficiently 
close to Mi (0) , u' a is also positive on (0, it) , and again u a intersects W\ at some point in (0, tt) . 
Let 

a 2 = sup{d G (wi(0), M p (0)) | For u\ (0) < a < a, u' a > on (0, tt) 
and u a intersects Wi at some point in (0,7r)}. 

Then a 2 < u p (0) since u p does not intersect Wi at all. 

By continuity, and since no u a is tangent to Wi, u 2 '■— u a2 intersects w± in (0, n] and u' 2 > 
in [0, 7r]. The first intersection of u 2 with w\ is either at tt or in (0, tt). In the first case u' 2 {tt) > 0, 
which implies {u 2 — wi)'( 7r ) > and so m 2 < w\ just before 7r, a contradiction. Hence «2 first 
intersects w\ in (0, 7r). Then by the definition of a 2 , it follows that u 2 = at some first t G (0, 7r]. 
We claim that t = n. If < t < n. then u 2 (t) = for otherwise, u 2 (t) < and this implies that 
u' 2 < just after t, a contradiction. So ^"(t) = — sint < which gives u' 2 < just before and 
after t, again a contradiction. This shows that u 2 (tt) = so that u 2 is a periodic solution which 
satisfies Wi(0) < ^2(0) < u p (0) and u 2 > in (0,7r). Further, u 2 intersects W\ in (0,7r), which 
implies that u 2 (tt) < wi(n). Let u^(t) = —u 2 (t + it). From Theorem |B| and Corollary ^, it again 
follows that Ui < u 2 < u 3 < m 4 < u 5 . m 

Remark 17 If the forcing term had less symmetry than cosine, but still was symmetric around 
it and 0, we could still construct the five periodic solutions by a similar shooting method. 

4.4 Verification of Condition [TT] for a range of parameters 



It is trivial to verify Condition |TT] for A > Ao and sufficiently small e. This is most easily seen 
by making the change of variables r = | , letting v (r) = u (t) , and considering the resulting 
equation 

v = v 3 — Xv + COS ST. (9) 
Then consider the phase plane for the limiting equation 

v = v 3 - Xv + 1. 

There are three equilibria in the v , v plane, p, = (ai,0) , i = 1,2,3, where a± < < a 2 < 0:3. 
The outer two, pi and P3 are saddle points while P2 is a center. There is a homoclinic 
orbit at p 3 , so that the left branch of the unstable manifold at p 3 is bounded, with v taking 
its minimum value at a point (a*,0) where ct\ < a* < a 2 , while the unstable manifold at p x is 
unbounded to the right. Suppose that v (0) = a G (ai, a*) . Then v > for t > as long as 
the solution is defined, and v — > 00 in finite time. Continuity with respect to e implies that 
for small enough e, the solution of (^]) with v (0) = a, v (0) = will increase monotonically and 
cross v = b before r = ^. This verifies Condition |11] for small e. 



Next we give an estimate on e for given A > Ao to ensure that Condition [TT] holds. This 
estimate of e is by no means best. 
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Lemma 18 For any A > X let b = 6(A) = y A + ^ and e\ = g|r, where T\ = 2\/b + yf\ for 
Ao < A < 4 and 



r A = 2v^ + ^ ln ^" 1 ^ 2(v / ^-v / ^ r ^ 

; £/ie solution 
there is a T G (0,T\) snc/i £/iat 



/or A > 4. Let v be the solution of ftQ) with v(0) = — x/\ and v(0) = 0. If < e < s\ in then 



Proof. Let 



v>0 in (0,T] C (0, ^-], and v(T) = b. (11) 

3e 



T = sup{r G (0, — )| v > 0,v < b in (0,r]}. 



Since v(0) = 1, it follows that T is well defined. Then for r G (0, T), cos(£t) > 1/2 and then (|9|) 
gives i; > f 3 — Af + 1/2. Multiply this inequality by 2?) and integrate over [0, r] to give 

(v) 2 > ^ 4 - At; 2 + i; - (^A 2 - A 2 - v^) = \{v 2 - A) 2 + (v + V\) (12) 
and hence, for r G (0, T], 



This implies that 



v> J~{v 2 -X) 2 + {v + V\). (13) 



, dv , , 

T< / - _ (14) 



and > (recall that v(0) = — n/A). 

i 



If A < A < 4, then T < -?=?dv = 2x/b+V\ = T x , which proves (|TT|) for < e < e x . 



Assume that A > 4 and write 

dv 



^ ^(v 2 - A) 2 + (v + y/X) 



h + h + h, 



where 



and 



h:= f - T iL= = 2(Vfc + VA-V / 2VA-2). 

It follows from (0) that T < T\ and therefore (|TTD follows from the definition of T and the 
assumption that < e < e\. ■ 
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Remark 19 

Condition 



Since lim 



2 a/2, it follows from Lemma [7^ that for sufficiently large A > ; 
77] holds for < e < v4/so, easy numerical estimates (assisted by a computer 



algebra program!) show that for any A > Aq, Condition 11 



holds ifO < e < \. 



5 Results for small e. 



5.1 Asymptotic form of periodic solutions as e —> 0. 

One of the key points made in the work of Angenent, Mallet-Paret, and Peletier is that as e — > 0, 
the solution which we have called u p , or u 3 , tends to the lower root, U (t) , of u" = in [0, |) 
and to the upper root, which we have denoted by U (t) , in (~, 7r]. (This was for the case A > A .) 
We will need this result here. Since one of the aims of this paper is to give elementary proofs of 
their results, with no reliance on infinite dimensional analysis or partial differential equations, we 
give a new proof. The result is obtained from the following lemma, where we make no assumption 
on A. For any A, the equation u 3 — Xu + cos (i) = has a smallest solution u = U_(t) which is 



continuous in any interval I k 
continuous in I k . 



(2fe-lW (2fc+l)7r 



with k an even integer. When k is odd, U is 



Lemma 20 For some even integer k let and be closed intervals with C int ( J/.) C Ik 
(where int denotes interior). Then for any 5 > 0, there is an Eq such that if < e < 6q, and 
if u is a solution of (1.1) with —b<u<0in Jk, then \u(t) — U_{t)\ < 5 in M k . If k is odd, a 
similar result holds, stating that solutions which are positive and bounded by b in J k are close to 
U in M k . 

If we recall that u p is an even function, we see that u p must be close to the lowest root U_(t) 
of u" = on any interval of the form [0, (1 — 5) |] as e —>■ 0. This is true also for any bounded 
solution with u (0) < u v (0) . 

Proof. Consider the case k even. The branch U_{t) is strictly negative in J k . Further, given 
a sufficiently small 5 > 0, depending on A and J k , the quantity 

/max = m & x f (t, u) 



is negative and the quantity 



/ min = min / (t, u) 

U(t)+5<u<0 



is positive. It follows Eq can be chosen so that if u(t) <U_ (t) — 5 for some t G M k and u' (t) < 0, 
then u crosses u = — b as t increases in J k) while if vl (t) > 0, then u crosses —b in J k as t 
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decreases. Similarly, if u(t) > U_ (t) + 5, then u will cross zero before t leaves Jj~, either forward 
or backward depending on the sign of u' . The case k odd is similar. This proves Lemma |20. ■ 



The remaining results in this section give more information about the asymptotic behavior 
of solutions as e — > 0. They will be needed in the next section. 



Lemma 21 Assume that X > An. Then 



(i) There are Mq > and Eq > such that if < e < eq and u is a solution of ([Jj 
satisfying u < — v/f over an interval [c, d], with — oo < c < d < oo, then u(t) < U_(t) + M e 2 in 

J £ := [c + j^e\ lne|, d — j^e\ lne|]. If u is a solution of (1.1) satisfying u < — y| m [c, d], taen 



n(t) > U (t) - M £ 2 in J £ . Here K = Ju(0) - y §, [c, d] := (— oo,d] i/c = -oo and d < oo, 
and other cases such as d = oo and c > — oo and both c = — oo and d = oo are defined similarly. 

(ii) If < e < Eq and u is a solution of (^) which is bounded in (—00,00) and satisfies 

A III [( . </ (17/;/ —..X. 2^ ( «• " 1_- -A- ■ "" " 

|n(t) - C/(t)| < M e 2 /or t G J e . (15) 
// instead u satisfies u > over [c, d] , then 

\u(t) - U(t)\ < M e 2 for t G J e . (16) 

Proof. It is easy to check that there are M > and Eq > such that if < e < Eq and 
U 2 = U + (M - 1)e 2 , then 

e 2 U 2 < U\ - \U 2 + cost 

for all t G (—00,00). Assume for contradiction that there are E G (0, £n) an d t G J e such 
that it > £/ + MqE 2 at the point t. Assume first that u'(t) > U!(t). Let w = u — 1/2- Since 
e 2 w" > (u 2 + uU 2 + U 2 - X)w in [c, d], w(t) > and n/(t) > 0, it follows that w > 0, «/ > and 
£ 2 w" > K 2 w for t G [t, d] . Then for £ < t < s < d we find that 



^ + y w 2 (s) - «; 2 (t) + -^( W ') 2 (t)j > j + ^n) 2 ( S ) - «; 2 (t) + ^{w') 2 {t) 

and so for t < t < T < d we obtain 



w(T) + \jw 2 {T) - w 2 (t) + ^(w') 2 (t) > [w(t) + ^w'Wef^. (17) 
Assume that d < 00. Evaluating ( |T7D at T = d and t — t and using d — i > -^e \ \ne\ gives 



w(d) + \/w 2 (d) - n; 2 (t) + ^- 2 (w') 2 (t) > Hi) + ^w'(t)}^ 
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Since wit) > e 2 and w'ii) > 0, it follows that the above inequality does not hold if e is small. 

Therefore, e can be chosen independent of t so that for < s < e , u'(t) < U!(t) if u(t) > 
U_+ M e 2 . Then let u{t) = u(—t) and Ui(t) = U^—t). Apply the same argument for it as above 
to get the same contradiction if e < s . 

If d = oo, then letting t — > oo, we see that the right side of ( |T7| ) goes to oo, contradicting the 
boundedness of the left side of ([171) . The other assertions of (i) can be proved similarly. 

(ii) Since U\ is the minimal solution of ([!]), it follows from (i) that to prove (|15|) it suffices to 
show that U\{t) > U — M e 2 in J e . We can choose Eq so that for < e < £o, 

e 2 U" > Ul - XUi + cost, 

where U\ — U_ — (M — l)e 2 . Assume that there is an e < Eo and t e J £ such that u\(t) < 
U_(t) — M e 2 . Then, as above, by considering (U\ — U\)(t) for increasing t if u[(t) < U[(t) and 
for decreasing t if Ui(t) > U[(i), we obtain that (Ui — u%)(t) — >• oo if t — >• ±oo respectively, a 
contradiction. This proves and (pT6|) can be proved similarly. ■ 



The next Lemma describes the asymptotic behavior of Ui, M5 and their first order derivatives 
as e — > 0. 

Lemma 22 Assume that A > Ao- M^e can choose Eq so small that if < e < Eq, then for all 
t G (—00, 00), 

\m(t) - U(t)\ < MqE 2 , and \u 5 (t) - U(t)\ < M e 2 (18) 

and 

Kit) - U'(t)\ < MqE, and \u' 5 (t) - U'(t)\ < M' Q e, (19) 
where M' Q is a constant independent of e. 

Proof. The inequalities in (|l^) follow from (ii) of the above lemma. To show ( |19|) we let 

z — ^2(^1 — JDi which satisfies 

e 2 z" = (3U 2 - X)z + 3e 2 Uz 2 + e a z 3 - Tf . (20) 

Further, z'(0) = and \z{t)\ < M for all t E (-00,00) and e G (0,e )- Multiply O by z' to 

get 

l -E 2 {z')\t) = l -{W 2 {t) - X)z 2 (t) - i(3f/ 2 (0) - A), 2 (0) -3^U(s)U:(s)z 2 (s)ds 

+ E 2 U{t)z 3 {t) - e 2 U{0)z 3 (0) - f z 3 (s)lf(s) ds + -E\z\t) - z\0)) 

Jo 4 

-U!'(t)z(s) + U! , (Q)z(0)+ [ z(s)U!"(s)ds (21) 

Jo 
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Since z, U_ and the derivatives of U_ are all bounded with the bounds independent of e, it 
follows that the right-hand side of fl2~T|) is bounded by a constant, say, |(Mq) 2 over [0, it], that is, 
e 2 {z'f(t) < (Mq) 2 for t G [0,tt] and s G (0,e ]. Hence holds for t G [0,tt]. Since u x and £/ 
are 27r-periodic and even functions, we see that (IB3) holds for all t G (—00, 00). ■ 



Lemma 23 Assume that A > Aq. TTien /or each small \l > t/iere zs an e M G (0, 1) such that if 



(22) 



< £ < and u is a solution with u> y ^ in [c, d] ; t/ien 



|ix — 1*5 1 + — — \u! — u' 5 \ < M%e ^ in [c + fi,d — fj]. 
2K 



Ifu<—\^ in [c, d], then 



u(t)-u 1 (t)\ + —\u'(t)-u' 1 (t)\ < M ie -T» m [c + fi,d-fi]. 



2K 



Here K is defined in Lemma [II] and Mi = 2(U(tt 



-) 

3) 



Proof, (i) Let w = M5 — u. Since w > and u > y ^ in [c, d], it follows that £ W > if iu in 

[c, d]. From Lemma |2^ we see that w < M e 2 in [c + /i, d — //]. We first assume that u>' > at 
to := c + /i. Then (|T7D holds for this if with t < t < T < d, and so for t < £ < ^ < d, 



w(t) + (1 - e"T V-t))l-\ w i(t)\ < 2w{T)e~TV-t) 



K 



(23) 



Since w(T) < U (tt) — y |, and e s M < | if £ is small, it follows that for t G [c — /j,, d — //], 

w(t) + ^\w'(t)\<M 1 e'^ T - t \ 



(24) 



Let T = d. Since d—t > fi for t G [to, d— //], (|22| ) follows immediately from ( pl"D for t G [c+/i, d— //]. 
Assume that w'(^o) < 0. Let t = sup{t G [to, d — /j] : w' < in [t , t)}. From what we just proved, 

we can assume that t = d — /i. Since u/ is decreasing in [c, d], io' < in [c — t], and w > in 
[c, d], it suffices to show that (p2|) holds at t — t . Integrating w"w' < over [c, to] gives 



«;(*„) + -l|u/(*o)| < [^(c) + ^/u; 2 (c)-w; 2 (to) + -^K) 2 ^o) 

which implies that ( p2]) holds at t = to- The proof of the inequality for tti is similar and therefore 
is omitted. ■ 

From Lemmas ETD, and E3] one immediately obtains a refinement of Lemma 121. 
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Corollary 24 Assume that X > Ao- For any small /x > ; i/xere is an > s?xca t/iat i/ 

< e < and u is a bounded solution of (fjj over (— oo, oo) satisfying u < ~\J\ over [ c >d], 
then for t £ [c + /x, d — /x] 

|u(t) -££(*) | +e|u / (t) - t/;'(t)| < (M + Mq + l)e 2 . (25) 
If instead, u satisfies u > y| over [c, d], then for t £ [c + /x, d — /x], 

|n(t) - C7(t)| + e\u'(t) - U'(t)\ < (M + M' + 1) £ 2 . (26) 



The next lemma is a refinement of Lemma 20 if A > Aq. 



Lemma 25 Assume that A > A . For any small /x > 0, there is an > s?xca i7ia£ i/ 
< £ < and u is a bounded solution of (Q/ 1 which satisfies u < — \/~X at some point in J kil : = 
j- (2fc-i)7r _|_£ (2fc+i)7r _£j some euen integer k, then $25] ) holds in M kpL := [ ^ 2fc ~ 1 - >7r +/x, ( 2fc + 1 ) 7r _^j i 

Similarly, $Zb\ ) holds in if u satisfies u > V~X at some point in J kfM for some odd integer k. 

Proof. To show the first part of lemma, from Corollary |24] it suffices to show that for e small, 

xx < ~\J\ m Jkfi- Assume that this is false. Then there is a sequence e n with lim^oo e n = and 

a sequence of t n £ J kfl with lim^oo t n = i for some i £ J kfl such that u n {t n ) = — y/\. Then by a 
phase plane argument we obtain that u n will reach b by time t + 0(e n ) as n — > oo, contradicting 
that u n is bounded. ■ 



Since u v is antisymmetric around 7r/2, from Lemma 25 we immediately get 



Theorem 26 Assume that X > X . For any small \x > t/xere is an e M > such that if 
< £ < e^, then ( p?3[ ) and ( flffi ) no/d /or u = u p in [— ~ + /x, | — /x] and [| + /x, ^ — /x] respectively. 

The next theorem is about the asymptotic behavior of n2 as e — > 0. 

Theorem 27 (%) Lei VL-j be the homoclinic solution ofV = V 3 — XV— 1 sxxc/x that lmv^-toc, V-\(t) - 
C/(7r) and V'-i(O) = 0. Then lim e ^ ^(tt) = VLi(O). For any given T > 0, n 2 (7r + er) — VLi(r) 
approaches zero uniformly for r £ [— T, 0] as e — > 0. (ii/ 1 For any smaK /x > there is an e M > 

snca £/xa£ if < e < e^, then $25] ) holds in [0, n — /x] /or n = n 2 . 

Proof. Let n(r) = n 2 (7r + er). Then n satisfies 

i3 = v 3 — Xv + cos(tt + er), v(0) = xx 2 (7r), v(0) = 0. 

Then (i) follows easily by a phase plane argument. To show (ii), we take T so large that 
V-i(—T) < — yf- Then from (i) it follows that if e is sufficiently small, then xx 2 (7r — eT) < 
Therefore, if e is small enough to satisfy eT < |, then, since n' 2 > in (0, ir), n 2 (t) < n 2 (7r — ^) < 

A for t £ [— |, 7T — |]. Hence (ii) follows from Corollary |4j. ■ 
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5.2 Further periodic solutions 



As pointed out earlier, the classical analysis of Duffing's equation, with results such as those 
]NM|, corresponded in a sense to taking A large. These solutions are bounded as A 



in 



oo. 



The solutions u±, 112, U3 found in section ||] all have minimum values (which are also their initial 
values) below — vA, while 114 and M5 have maxima above a/A. In this section we begin by finding 
solutions with oscillations, distributed reasonably evenly in (0, 2ir) , and with initial values in the 

interval ^— j . This is for a fixed A and small e. Then we consider the behavior of those 

solutions if we keep the number of oscillations fixed (together with A) and let e —* 0. It is 
found that the oscillations collect near odd multiples of |, forming internal "layers" between U 
and U_. (See Figure 6.) The result is related to one in [ |HM2|| about a similar nonlinear forced 
oscillator (derived in [pOJ|| ), where the nonlinearity was quadratic rather than cubic. However 
the Duffing equation has a richer collection of solutions than the equation in [[H1V12 . 




A question of interest is to what extent the results in this paper depend on the symmetry of 
the cosine function. This is addressed further in section 0, but the techniques in the Theorems 



28 and 31 of the current section do not use the symmetry of cosine around ~. These should apply 



to more general forcing functions, as we will discuss in future work. On the other hand, the proof 
of Theorem ^ below shows how use of the symmetry around ~ can greatly simplify some proofs. 
The key fact is that if u is a solution of (D) such that u (|) =0, then u (| — t) = —u (| + 1) for 
all t. 
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In some of the following results we will switch back and forth between the original scaling 
and that in the slow equation (^) , with t — er, u it) — v (r) . For any k G (0, 1] , the equation 



V = V s - XV + k (27) 

has a unique solution, forming the homoclinic orbit, with V (±oo) = U (arccos k) and V (0) = 0. 
We denote this solution by V K . Let 

A = inf{A > A : Vi(0) < U (tt)} 

Thus, for A > A, the homoclinic orbit at r = (k = 1) has a minimum below Uq (tt) . It is easy 
to show that A < 3. 



Theorem 28 Suppose that A > A. Then for any N > there is an eo such that if < e < 6q, 
and if M is any positive integer with M < N, then ([!]) has a periodic solution with exactly M 
local maxima in (0, n) , the last being at it. Also, these periodic solutions have a local minimum 
at t = and satisfy —\/~X < u(0) < U (n) and U (ir) < u(ir) < U(n). 

Proof. We will need the following lemma, which is similar to Theorem 2 in [rlM2|. 

Lemma 29 Let u be a solution of ([][) with u'(0) = 0. If ti and t 2 are two successive minima 
(maxima) of u with < ti < t 2 < n, then u(t 2 ) < u(ti). 

Proof. We only show the case that u(ti) and w^) are successive maxima of u and the 
other case can be proved similarly. Assume that u reaches its maximum between t\ and t 2 at 
t = t m . Then multiply ([[]) by u' and integrate over [t_,i + ], where t\ < £_ < t m < t + < t 2 and 
u(t_) = u(t + ), to yield 

-e 2 [{u') 2 (t + ) - (u') 2 (t_)]= u'(t)costdt= ( / +/ )u'{t)costdt 

v t v t J 

u(tm) 

(cost-(-u) — cost + (-u)) du > 0. 

«(*-) 

Here we use that cost is decreasing in (0,7r), and t_(n) and t + (u) are the inverse function of 
u = u(t) for t G [ti,t m ] and t G [i m , ^2] respectively. It then follows that u'(t-) = before 
u'(t + ) = 0, which implies that u(ti) > u(t 2 ) as required. 
■ 

For A > A , recall that U(t) < U (t) < U(t) are the three solutions of u 3 (t) — \u(t) +cost = 0. 
For given A > A, let 5 = Uo(tt) — Vi(0). We also need the fact, easily proved, that if V K is the 
homoclinic solution of (|5]) defined earlier, for < k < 1, then the minimum values V K (0) are 
decreasing with respect to k. 

Also, for each a, let v a be the solution of (|) such that v (0) = a, v' (0) = 0. Let a = 
U (tt) - f . 

Choose > 0. When e — 0, v& is periodic, and so has A^ + 2 local maxima in some interval 
[0, T/v] . For sufficiently small e, V& still has at least A^ + 1 local maxima in (0, ^) . Hence U& 
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has at least N + 1 maxima in (0, tt) . Suppose that these are at < ti < t 2 < ... < t^+i < tt. 



Further, from Lemma |29| we see that (t i+ i) < (U) for < i < N. 



If any of the maxima ti for i < N are such that Ua (ti) < U_ (ti) , then because U_ is increasing 
in (0, tt) , we must have u < U_ in (ti, tt) (if the solution exists out to tt), and so could not have 
any more maxima. This shows that for 1 < i < N, the maxima of must lie in the interval 
\Uq (t) , U (£)] . However, if u' = u" = at some t in (0, tt) , then u'" = — sin (t) < 0, so this 
could not be a maximum. Hence the first N maxima of must lie in (U (t) , U (t)) . 

The solution must have a minimum at so = and further minima s\, ....sn £ (0, tt). Also 
a > u (s\) > u (s 2 ) > ... > u (sjv), by Lemma Therefore, u& (U) > + | > (s^ + | 

fori = l,...,iV. 

Now decrease a from a. Since u" < at a maximum, the maxima are continuous in a, and 
so remain above Uq as long as they exist. If a = —b, u a has no local minima. But as long as 
there are, say, M < N local maxima in (0, rr) with u > Uq and with a < a these maxima and 
their intervening minima are separated, in that u a (ti (a)) — u a (s.j (a)) > | for 1 < i < M. This 
means that the number of maxima cannot decrease until one crosses t = tt. This must happen 



successively for each maximum, which proves Theorem 28 



Remark 30 With further estimates it can be shown that for small e there are 2tt periodic solu- 
tions with as many as y maxima in (0,7r) ; where K > is independent of e. 

In the next theorem we shall describe the asymptotic behavior, as e — » 0, of 27r-periodic 
solutions of ([l]) with m maxima in (0, tt] and with u (0) < 0, u (ti) > 0. The result, informally, 
is that if m is kept fixed as e — ► then the internal maxima and minima all tend to |, while 
near and tt the solutions tend to homoclinic orbits of the appropriate limiting equation. Thus 
these solutions have "spikes" at and tt, and internal layers where they are close to appropriate 



heteroclinic orbits, near ^. In theorem we obtain 27r-periodic solutions without spikes at 
and tt, using the symmetry of the cosine. In forthcoming work we expect to consider other 
multilayer solutions, combining spike and non-spike behaviors at multiples of tt, without reliance 
on symmetry. 



Theorem 31 Suppose that A > A. For an integer N > 2 choose Eo as in Theorem [|<|. Let m 
be an integer with 1 < m < N, and for each e G (0, Eo) let u = u £ be a 2tt periodic solution of ([|) 
satisfying u'(0) = u'(tt) = and —y/X < u(0) < Uo(0), and having exactly m maxima in (0, tt]. 
Let t\, ■ ■ ■ ,t m and s\, • • • , s m _i be the maxima and the minima of u in (0, tt] respectively such 
that < t 1 < si < t 2 < s 2 < ■ ■ ■ < t m _i < s m _! < t m = tt. Let V\ be the homoclinic solution 
of v = v 3 — Xv + 1 such that lim T ^± OQ v(T) = U(0) and v(0) = 0. Suppose that 5 is a positive 

number with 5 < U(0) — Let V-i be the homoclinic solution of v = v 3 — Xv — 1 such that 

lim^-too v(r) = U_(tt) and i)(0) = 0. Let Vq+ be the heteroclinic solution of v = v 3 — Xv with 
lim^oo v(t) = [7(f), lim^oo v(t) = C/(|) and v(0) = [7(f) - 5. Also, let V - (t) = -V + (-t) . 
Then, as e —>■ 0, 

(i)u £ (0)^V 1 (0),U £ (7T)^V^ 1 (0). 
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(%%) for each 1 < j < m — 1, tj — > \, Sj — > | ; u e (tj) — > £/(f ), and — ► C/(f)- 

(Hi) for any given T > 0, \u e {er) — Vi(r)| + \eu' £ (er) — Vi(r)|) — > uniformly for r G [0, T], 
and \u £ (tc + er) — V-i(r)\ + \eu' £ (n + er) — V-\(t)\) = uniformly for r G [— T, 0]. 
Also, for 1 < j < m — 1, let 

Tj = sup{t G (tj, 7r) : m £ > m 5 — 5 in (tj, t}} 



Sj = supjt G (sj, tt) : u e < U\ + 5 in (sj, t]}. 



Then < Tj — tj < ^|lne| and < Sj — Sj < l§|lne| 7 where K is defined in Lemma 21 



Further, for any given T > 0, as e — > 0, |« e (3} + er) — V + (r) | + \eu' e (Tj + er) — Vb+ (r)| ^ and 
1^(5^ + er) + Vb- (r) | + \eu' e (Tj + er) + V - (r) | -> 0, uniformly for r G [0, T] . 



Finally, for any small fi > t/jere zs an e M > stzc/i £/ia£ if < e < e^, then $ZB) and 
hold in {pi, f — /jl] and [~ + fi, tt — /jl] respectively. 



Proof. 

In the proof we shall suppress the dependence of u on e. We first show that lim £ ^ M (0) = 
Vi(0). Assume that this is false. Since w(0) is bounded, there is a sequence e n with lim^oo e n = 
such that as n — > oo, the corresponding u(0) approaches (3$ G [— V% £/o(0)]. For simplicity we 
assume that lim e ^ ^(0) = Po- Let z(r) = u(er). Then z satisfies 

z = z 3 - Xz + coser, z(0) = u(0), i(0) = 0. 

Let z be the solution of 

i = - A^ + 1, z (0) = p , z (t) = 0. 

Since lim £ ^ z(0) — p , it follows from the continuity of solutions with respect to parameters that 
for any given T > if Zq is defined on [0, T], then lim e ^ Z ( T ) — z o( T ) uniformly for r G [0, T\. 

Assume that p < Vi(0). Then there is a T > such that Zq(t) > for r G (0, T ) and 
lim T ^ T - Zq(t) = oo. It follows that for e sufficiently small z crosses b, which is impossible. 

Hence P\ > Vi(0). Then z is a periodic function with a period T > and so z has m + 1 
maximum in the interval (0, (m + 2)To). Hence by continuity, for e sufficiently small, z(r) also 
has m+1 maximum in the interval (0, (m + 2) To) , which implies that w has m+1 maxima in 
(0, e(m + 2)T ) C (0, tt] for e sufficiently small, contradicting the assumption on u. Therefore, 
lim £ ^ ^(0) = Vi(0) and so (hi) and the first part of (i) follow. The rest of (i) and (hi) can be 
proved similarly. 

We next show that lim e ^ ii = f • Suppose not. Since t\ G (0, tt], there is a sequence e n with 
lim n ^oo e n = such that lim^oo t% — i ^ | for some i G [0, tt). Again, for simplicity, we assume 
that lime^o t\ = t. 
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Since u{tt) — > VLi(O) as e — > 0, it follows from Lemma ^ that t ^ tt. We suppose now that 
< £ < 7T. Let z(t) = u(tx + er). Then z satisfies 



z = z 3 - Az + cos(ti +er), z(0) = u(t 1 ), z(0) = 0. 



2A 
u(ti) 



For 6 = 6 (A) = \/ \ + ^ (as in Lemma 0), let 

dy 



and 



6 ^(Z/ 4 - - KV 2 - «(^i) 2 ) + % - w(ti)) costx ' 



T = sup{r G (0,T 6 + 1) : z > -b, z<0 in(0,r)}. 



Observe that U (t\) > z(0) = u(ti) > U (ti)—S if e is sufficiently small. Hence, T& is defined and 
Tt < oo; z(0) > and so T is well defined. Assume that e is so small that ti+e(Tb+l) < tt. Then 
on (0,T) we have z < 2 3 -Az + costi and so {zf > \(z A - z A (Q))-\(z 2 - z 2 (0)) + 2(z- z(0)) cos^. 
We see that T < and i(T) < 0. Therefore, by the definition of T, it follows that u(ti + eT) = 
z(T) = —6, which is impossible since u{r) > —b for all r. 

We assume now that i G [0, |). Note that u(ti) < U[t\). We shall show that lim £ _».o u{ti) = 
U(t). For if this is false, then the boundedness of u(ti) implies that there is a sequence e n , with 
e n — > as n — > oo, such that lim^oo u(ti) < U(t). It then follows by arguments similar to those 
above that u has more than m maxima in (0, tt] for large n, which is a contradiction. 

We now claim that \T\ — 1\\ < %| lne| for e > sufficiently small, where T x is defined in the 
statement of this Theorem (31). Assume that this is false. Let w = — u. Then in (ti,Ti), 
w > and e 2 w" > i^ 2 w. Since u' < and u' 5 > just to the right of t 1) we see that w' > and 
w"w' > (K 2 /e 2 )ww' in (ti,Ti]. Then, as in fll7|) (where w was slightly different), we again get 



l r(T)+ x j w ^{T)-w 2 {t) + ^{w') 2 {t) > [w^ + ^w'^etV-V. (21 



We 



assume now that t > 0. Since lim e _^ U'(ti) = U'(t) = 3 fj2^_ x , we use C PL9[) to obtain 

w '(t x ) = u'Atx) > U'itt) - M'e > -U'(t) - M'e > , - ' r > 0. 

ku 5V v - v i; o _ 2 w o 4(3f/ 2 (t) - A) 



This estimate gives a lower bound on w' (ti) as e tends to zero. Evaluating (p8|) at t = t x and 
T = t x + ^| lne| and letting e — > shows that the right side of (fffi ) goes to infinity, while the 
left side is bounded. This contradiction proves the claim if £ G (0, |) . 
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We now assume that t — 0. We observe that 

w% + -|| lne|) > u' B (ti + -|| lne|) > + -|| lne|) - M' s 

sin(^|lne|) gljnfj 
4(3£/ 2 (0)-A) 2vr(3f/ 2 (0)-A)' 



Then evaluating (p8|) at i = t% + jA lne| and T = t± + M \ In e| and letting e — > 0, the same 
contradiction will be obtained. This shows that \T\ — t\\ < |§| lne|, as claimed. In particular, 
T\ — t\ — > and so T\ — > i. 

Continuing with our proof that t = ~, and under the assumption that this is false and 
t e [0, |), let z(t) = u(T\ + er). Then z satisfies 

2 = 2 3 -Az + cos(T 1 + £r), z(0) = u{Ti) = u(Tx) — 8, z{0) = eu'(T 1 ). (29) 

We shall show that lim £ ^o^(0) = io(0), where zq is the homoclinic solution of 

z = Zq — Xzq + cos(F) (30) 

with zq(0) = U(t) - 5 and i (0) < 0. 

Suppose that lim £ ^ z(0) ^ io(0). Since i(0) < and i(0) is bounded, which is easily verified, 
there is sequence of values of e, and corresponding solutions u with corresponding z, such that 
z (0) approaches a number a ^ io(0) with a < 0. We first assume that a > zq(0). Let z\ be the 
periodic solution, with period T, to 

z x = z \ - Xzi + cos(t), z 1 (0) = U(f)-6, z 1 {0) = a. (31) 

Since Ti — > t, 2(7") approaches 2i(t) as e — > uniformly in compact intervals of r, and it follows 
that for sufficient small e, z oscillates more than m + 1 times in [0, (ra + 2)T], and so u has more 
than m maxima in [0, tt], a contradiction. 



Hence we can assume that a < Zo(0) and again assume that Z\ solves (pi]). Then there is 
a T > such that ii(r) < for r e [T, 0] and z\{T) > b, which implies that for e sufficiently 
small, u(t\ + eT) > b, again a contradiction. 

Therefore, as e — > 0, i(0) — > io(0) and z goes to z uniformly in any compact interval. 
Since zq is homoclinic to U(t), it follows from continuity that after Ti, u will return to the any 
given neighborhood of U(t) before t — T\ + Me for some M independent of e. Then, since 
m = 2, u increases and stays close to U till tc, which implies that lim e ^ M ( 7r ) = U(ir) and so 
u{t2) = u(ir) > u(ti) for e small, contradicting Lemma |29|. This shows that i = | for m = 2. 

Then, zo is the heteroclinic solution connecting £/(§ ) as t — > — oo to ££(§) as £ — ■> oo. So, 
as e — >■ the point (£,«(£)) reaches a point as close to (|, C/ f|)) as we like, and then Lemma 



29| implies that lim e ^ s i = f an d lim e _» w (si) = ££(!)• For t e [ti, 7r], w (t) < u (ti), again by 



Lemma 123. Further, u(t) remains close to U_(t) until t is close to ir, for otherwise there would 
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be similar contradictions to those obtained above. Hence, u must, after its last minimum, follow 
U_ until close to it and then , since its last maximum is above Uq(tt), follow V-\. The bound on 
Sj — Sj follows in the same way as the bound for Tj — tj. For the final statement in the Theorem, 
we use (i)-(iii) and Corollary E3. This proves the theorem for m = 2. The proof for m > 2 is 
similar. ■ 



The solutions discussed in Theorem |31| have upward and downward pointing "spikes" at each 
multiple of 7r, corresponding to homoclinic orbits in the phase plane, as well as "layers" at odd 
multiples of | corresponding to heteroclinic orbits. In the next theorem we use the symmetry 
of cosine around | to give a quick proof that there are also solutions with layers but without 
the spikes. In our future work we expect to show that these solutions exist without reliance on 
symmetry, but it is also valuable, we believe, to show how quickly a proof can be obtained in the 
symmetric case. We need a preliminary lemma. 

Lemma 32 Assume that A > Ao- Let u be a solution of ([|) with u(|) = and u'(^) = (3. 

(i) If (3 < then either u' < in [| ,ir) or there is t Q £ (f,7r) such that u(t ) = U_(t ) 
andu'{t) < in [f,io]- 

(ii) Let H(u) = \u 2 — u 4 /2. Assume that (3 < — V^pM l _ jff 1 g (|,7r) is the first time that 
u' = (if such a point exists), then u(ti) < U_{^)- 

Proof. From we have 

e 2 {u')\t) =e 2 {u')\-)-H(u) + 2 f v! cos sds. (32) 
2 J f 

(i) Let t = sup{t E (| ,7r) : v! < and u > -y/X in }. It follows from (||D that for 
t G (| , to), 

e 2 {u') 2 {t) > e 2 f3 2 - H{u) > e 2 (3 2 - y > 0, 

from which the assertion in (i) follows. 

(ii) From (0), we first have — \/~X < u(ti) < 0. Evaluate fl3"2| ) &t t = t± to give 

rh 

-2 r>2 



= e 2 (3 2 -H(u (h)) + 2 u' cos (s) ds 
so that 

H(u (ti)) > e 2 l3 2 . 
This implies that u (ti) < U_(n), as desired. ■ 
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Theorem 33 For any given integer m, there is an e m > such that for e G (0, e m ), has a 
solution u such that u'(0) = u'(tt) = 0, u(^) = 0, u'(|) < 0, and u has m minima and m maxima 
in [^, 7r]. If we denote these minima and maxima by Sj and tj respectively, with 1 < j < m, then 

TT 

- < S 1 < t X < S 2 < t 2 < ■ ■ ■ < S m < t m = 7T, 

71 71 

limsi = • • • = lims m = -, limu^sj) = • ■ ■ = limu(s m ) = U(-), 

e^O e— »0 2 e^O 2 

TX - 7T 

limti = • • • = limi m _! = -, limw(ti) = ■■■ = limu(t TO _i) = U(-), 

e^O e— >0 2 £^0 e^O 2 

u'(t) > 0, C/(|) < u(t) < U(n) for s m <t<n. 

Proof. Let up denote the solution u of ([[]) satisfying w(f ) = and u'(tt) = [3. Choose 

^ l e ( — "^e' ~ y / H ^ n —) and /?2 < — By a phase plane argument, there is an e m such that 
for e < e m , the solution up x has at least m minima and m maxima in (|, 7r) with all its first m 
minima lying between {7 and C/ an d all its first m maxima lying between Uq and U . Denote 
the m-th minimum and maximum by s m and t m . We consider the change in these minima as (3 



decreases from (3\. Part (ii) of lemma 32 shows that all the minima in (|,7r] lie below the line 
u = U_{tt), and so they can neither pass through t = tt nor disappear at the middle branch Uq of 
u" = 0. Since up 2 does not have any minimum at all, it follows that as we decrease f3\ to /?2, all 
the minima will disappear by crossing the lower branch of u" = before t — tt. 

Let = inf {/i|s m is defined continuously for (3 G (/3,Pi) as the m th minimum after |} . Then 
u'g 3 (s m ) = and U/3 3 (s m ) = U_(s m ) . For (3 — /5 3 > small, the m-th maximum t m {(3) of up exists. 
Clearly lim^^ t m (fi) = s m (fc). As before, if u p (t m ) = U(t m ), we get u'^t) < for all t - t m 
small, contradicting the definition of t m . Hence up{t m {(3)) < UJt m ((3)). Therefore, as we raise 
(3 from /?3, t m (/3) has to move toward to tt under the lower branch U_ of u" — 0. This maximum 
tm (P) cannot disappear by merger with another minimum in [f , 7rl , since this minimum would 
have to lie above up(s m ), contradicting Lemma Since up^tm) > Uo(tt) > U_(tt), there 

must be a point of discontinuity of the mth maximum before reaching /3i, and so there is a 
{3\ G (fa, Pi) such that t m (P) is continuous in [/3 3 , /? 4 ], t m {P^) = n and M^ 4 (t m (/?4)) = 0. Then 
up 4 is the desired solution. The asymptotic formulas stated in this theorem can be proved in a 



similar way to that of Theorem 31 



Figure 7 shows graphs of two of the new solutions, with one or three "layers" near |, for 
e = 0.5, A = 3, together with m 2 and the function uq (t) . It is of interest that these solutions 
behave in a very different way from the solution u p discussed earlier. While u p increases near | 



and decreases near ^r, these solutions behave in the opposite way. 



The existence of these solutions, and many more, with different numbers of oscillations 
between different zeros of cosx, was conjectured independently by H. Matano (private commu- 
nication). This was for a more general class of equations, but the result was limited to a finite 
interval, so that "chaos" was not involved. Professor Matano calls the solutions with the op- 
posite behavior from u p "up wind" solutions, a graphic term which we adopted in the caption 
of Figure 8. Similar solutions were obtained for a different equation by Nakashima [ JNAK|| . We 
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wish to thank Professor Matano for helpful correspondence, and in particular for sharing his con- 
jecture with us. While the results in this paper, about periodic solutions with multiple internal 
layers, were obtained independently, our thoughts concerning non-periodic solutions with more 
than three internal layers near odd multiples of | were previously somewhat vague, and we have 
been inspired to pursue this topic further by Professor Matano's conjecture. We expect that a 
proof of this conjecture using shooting methods, at least for some class of forcing functions, will 
be part of our next paper. The methods proposed by Matano for obtaining these solutions are 
very different from ours. 



Nakashima, in IINAK 



has also studied the stability of the oscillating solutions for her equa- 
tion, including the dimensions of the unstable manifolds (Morse Index). This has implications 
for the dimension of the global attractor for the problem (0) — (H) . Results of this type were also 



obtained in [jEM2], for a different equation, though they were not stated in these terms. We 



expect to study these topics for the current problem in future work, using the methods of [[HM2 
which are very different from the approach taken by Nakashima. 



Figure 7 

Graphs of two "up wind" solutions, together with w 2 and Uq. 

We end this section with a numerically computed solution suggesting how the different kinds 
of solutions can be combined. In this figure, e = 1, which explains why the oscillations of the 
new solution are so small. Results about this kind of solution will appear in part II. 
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Figure 8. 



Up and a new solution with three "up wind" sections 



5.3 Isolation and stability of u 



In this section we show how a shooting method can give results similar to those in 
without the use of infinite dimensional analysis or abstract dynamical systems. Using standard 
ode methods we prove the linearized stability of the three solutions found in section ^| with 
respect to (H) — (Q), and the isolation of u p from other solutions satisfying u' (0) = u' (kn) = 0. 



Full nonlinear stability follows by standard methods, laid out explicitly for this problem in fBF 
However, stability is inherently an infinite dimensional problem, and we certainly do not claim 
to prove (nonlinear) stability without the use of functional analysis. 

Theorem 34 Suppose that A > Ao- Choose r\ with < r\ < \. Let 

>-9> 



[i = / sin s ds. 

J (l—ri)n 

Choose 5 > so that 

(0 6<U(£*)-U{0) 
(ii) If U (t) — 5 < u (t) < U (t) for some t, then \u 3 — Am + cost| < 

Let Ss denote the set of points (t,u) such that < t < tt, —b<u<b, and 

-b < u < U_{t) + 5 if < t < ri7T 

U (t) - 5 < U < b if (1 - n) 7T < t < 7T. 

Then for sufficiently small e, u p is the only 2'K-periodic solution of ([I]) with u' (0) = u' (n) = 
whose graph over [0,7r] lies in Ss- Further, M p |[o i7 r] is a stable attr actor for the problem (H) — 
(^|) .with L = 7r 

Also, suppose that Ss is extended to [0, 2tt] by the operation (t, u) — > (2n — t,u) , and then 
periodically to [0, oo). Then for any positive integer k, u v is the unique solution of (P satisfying 
u' (0) = u' {kn) = whose graph lies in the extended set Ss, and its restriction to [0, kn] is a 
stable attr actor for @ — ([D with L = kn. 



To compare this with the result in [|AMPP|| we notice that the width of the vertical strip in 



Ss where the solution can increase from near U_ to near U is (1 — 2r x ) 7r, and this is free to be 
chosen within the constraint < r\ < |, By contrast, ||AMPP|| state only that there is some 



strip, with width independent of e, which contains the "internal layer" of the solution we have 
denoted by w P) and no other solution satisfying the boundary conditions has a jump upward 
within this strip. 
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The proof in [[AMPP|] is by a detailed construction of u p using sub-and super-solutions. 
It also uses abstract results from dynamical systems (to get the uniqueness.) Our proof of 
uniqueness is more direct, and starts with the existence of u p as given in Theorem |6|. 

Proof. To prove Theorem ^ we consider the variational equation and initial conditions 
satisfied by v = These are 

eV = (3 M 2 - A) v (33) 
v (0) = 1, v (0) = 0. 

We will also be concerned with w = u', which satisfies 

e V = (3m 2 - A) w - sin t (34) 
w (0) = 0, w' (0) = u" (0) 

We observe that w' (0) > when u (0) G (U (0) , 0) . Multiplying (|3|) by w and (g|) by v, 
subtracting, integrating by parts and using the initial conditions on v and w, we obtain 

r* v (s) 

wv ' - vw'L = - u " (0) + / sin s ds. (35) 

Jo £ 



Lemma 35 If u is a solution with u'(0) = which remains in Ss on [0, ir] , then v > on [0, tt] 
and v' (tt) > 0. 

Proof. Because u < U_(t) + S over the interval [0,ri7r], and therefore 3m 2 — A > in 
this interval, (|33]) implies that v grows exponentially large. More precisely, there are positive 
numbers K\ > 1 and 7, independent of e, such that v > K\e^ 1 in [0, tyk\ . Also, u" (0) = O (\j 
as e —>■ 0. Therefore, for small e, the right side of (^) is positive as long after t = r\-n as v is 
positive (up to t = n). 

We now show that v! ((l - ^) n) > 0. If not, then u' < on ((l — y) tt, rr) , since u" < 
when U — 5 < u < U, and U' > in (|,7r). Therefore m (7r) < u ((l — ^) 7r) . From (i) it 
follows that u (n) < U (tt) — 5, contradicting the assumption that u remains in Ss- Hence, 
u' ((l — y) 7r) > 0. This implies that v! > on (0, (l — y) tt], for if not, then u would have a 



minimum in this interval, and this minimum would lie above u (0) , contradicting Lemma 29 



Now suppose in (|35D that v = somewhere in (0, (l — y) tt] . Then the right side of (^) is 
positive, while the left side is negative. Hence, v > on [0, (l — y) 7r] . 

It is possible that v' becomes negative somewhere in [0, tt] . Indeed, numerically this is seen 
to happen. However, we will show that v' ((l — y) 7r) > 0. Suppose that v' ((l — y) 7r) < 0. 
Then v'((l — ri)7r) < 0, for otherwise v' would be positive and increasing on [(1 — r\)ir, (1 — y)vr), 
because (3m 2 — A) > there. Thus v is positive but decreasing on [(1 — ri) tt, (l — y) 71"] . From 
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(]3~5|) we obtain 



''(( 1 -t)')^(( 1 -t)')-(( 1 i -^)')-'(( 1 -t)') 

= — w" (0) + / — — sin s (is + / 2 sin s c?s 

/. n , \ /'( 1_ ^) 7r sin s , 
>«((l-y)T) / "71-^- 

v 2 / y (1 _ n)7r 

Since v! ((l — y) 7r) > 0, we get a contradiction from condition (n) in the statement of the 
theorem, and so v' could not remain negative on [(1 — r{) tt, (l — y) 7r] . Therefore, at (l — ^) 7r 
we have both v and u' > 0, and since 3m 2 — A > on [(1 — y) 7r, 7r] , both remain positive out 
to tt. This proves the lemma. ■ 



In particular, this applies to the solution u p . It follows by standard stability arguments [ |BF| 
that the solution u p is a stable attractor for the problem (fj) — (^) on [0, tt]. 

To prove that u p is the only 27r-periodic solution which remains in Ss, it is convenient to 
truncate the nonlinearity in by replacing u 3 — Am by b 3 — Xb for all u > b, and similarly, if 
u < —6, replace u 3 — Xu in ([!]) with Xb — b 3 . This means that all solutions u a exist on [0, n] , 
and we can consider u' a (ti) to be defined continuously for all negative a. We also note that 
if u a is a solution which does leave the region [—&,&] , then from the point where \u a \ = b, \u' a \ 
continuous to increase and cannot satisfy u' a = 0. 

Recall that there is a unique cti such that u ai is 27r-periodic and lies entirely below 
this is the "minimal" bounded solution, and its graph does not lie in S$- In considering the 
possibility of a second 27r-periodic solution, besides u p , which lies in S$, we need only consider 
a G (ai,U_ (0) + 5) . 



The solution u p remains in the interior of S$ on [0, it] , and the same is true for u a if \ct 



Otr. 



is sufficiently small. However as we raise or lower a from a p , we reach values where u a leaves 



Ss in [0, 7r] . Let / 



(3, a 



be the maximal interval containing a p such that u a remains in S$ on 



[0, 7r] for all a E I. I is well-defined and closed because Ss is a closed set. 



Lemma 36 If a ^ /, then u a does not remain in Ss on [0,7r] . 

Proof. : We claim that leaves Ss at (n, U (n) — 5) and exits Ss at (it, b) . If not, 
then one of these solutions is tangent to the boundary of Ss at some t < tt. For example, a 
tangency could occur at ((1 — r%) tt, U ((1 — r{) tt)) . But then, no matter what the slope of u is 
at this point, a phase plane argument shows that for sufficiently small e, u must cross —b on 
one direction or the other, and so nearby solutions also leave Ss before t = tt, contradicting the 
definition of a or (3. Similar considerations apply at any other possible tangent point in [0, tt). 
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Consider the case of us. By Lemma |35], -§^u a (tt) | a _g > 0. Hence, for a < (3 and close to f3, 



u a (k) < U (it) — 5, and we cannot, as we lower a, find a lower a where u a remains in Sg and 
u a (tt) = U (it) — 5, for at the first such point we would have -§^u a (it) < 0, a contradiction to 



Lemma 35. Similar remarks apply to u&, completing the proof of Lemma 36 



The proof that if u = u p , then v > on [0, ir] and v' (it) > 0, shows that in some neighborhood 
of dp, (a — a p ) u' a (n) > 0. Now, if there is an a G (a p , a) with u' & (tt) = 0, then choose the 

smallest such a. By Lemma [J5] we get -j^u' a (it) > at a p and a, a contradiction because these 
are adjacent zeros of u' a (7i). This contradiction can be reached similarly in the interval (j3,ap^, 
completing the proof of the uniqueness of u p among solutions with period 2tt which remain in 

s s . 

To extend the uniqueness and stability statements to larger intervals we note that starting 
with v (tc) > 0, v' (n) > 0, the same analysis allows us to show inductively that v remains 
positive, and v' is positive at any multiple of it. The arguments about uniqueness and stability 
can then be extended to [0, kn] completing the proof of Theorem B3. ■ 



5.4 Sensitivity with respect to initial conditions 

As this paper is already quite long, we will content ourselves with a few remarks. Up until 
now there has been no mention of "horseshoes" in this paper, or of Poincare maps, because our 
technique is to follow complete solutions of the ode, rather than to take snapshots at regular 
intervals. The results, however, are related to standard dynamical systems concepts such as 
horseshoes and sensitivity to initial conditions. 

In Theorems ^ and [15] we obtain a weak kind of sensitivity to initial conditions. For each 
sequence there is a corresponding solution, but the relation established is not 1:1. There could be 
an interval of initial values a in which u a intersects the same sequence of w^. This corresponds 
to a so-called "topological" horseshoe, without the hyperbolicity that was a key feature of Smale's 
original derivation. (See | |G.H]| for general discussion and references.) It is noted, however, that 



the solutions in Theorem [12| corresponding to the sequences of all odd integers and of all even 
integers, are our periodic solutions u\ and u 5 , and for these, uniqueness of the correspondence 
is established. The difficulty is the stability, or at least hyperbolicity, of solutions «2, «3, and 



In Section |5.3| , hyperbolicity was established for m 3 when e is sufficiently small. Theorem |33 
implies that any solution which corresponds to a sequence with no 2 or 4 is isolated from any 
other such solution. Therefore, in Theorem if we consider a sequence with no 2 or 4, then the 
infinite intersection of closed intervals used in the construction contains exactly one point. This 
means that the set of solutions found in Theorem [12] corresponding to sequences chosen from 
the set {1, 3, 5} is, for sufficiently small e, in 1:1 correspondence with the set of allowed symbol 
sequences. (The rule that 1 and 5 cannot follow each other must still be obeyed.) Hence the 
desired degree of sensitivity to initial conditions, in which any small perturbation of the initial 
condition leads to a deviation from the given sequence, is achieved. 
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In the paper [ HM2|| , a similar result was obtained. There, however, the analysis was not 
only for stable solutions. Indeed, only one of all the solutions found in that paper is stable in 
the linearized sense. However they are all hyperbolic. We believe that a similar analysis will 
allow us to study linearizations around solutions U2 and and prove their isolation as well. 
However we will not attempt this analysis here. 



5.5 Bifurcation in A 



We saw in the Introduction that at some A& G (0, Ao] new periodic solutions appear. While we 
know that for A > Ao there are at least five solutions with period 2ir, we have not discussed the 
nature of the bifurcation. It is possible that the number of solutions goes from one to three at 
\b and later increases to five. Stability considerations do not rule this out. 

One way of studying this numerically is to consider the graph of G (a) = u' a (tc), since we 
are only considering periodic solutions with the properties u' (0) = u' (tt) = 0. If the bifurcation 
is "pitchfork" , then the graph of G will qualitatively resemble that of the function a 3 — fia as 
the parameter /i changes from negative to positive. (The zeros of G will not be at a = 0.) 
However a numerical study of the function G quickly suggests that this is not the case. In 
Figure 9 we show the resulting graph of G for e = 1 close to the bifurcation point, A = 1.023. 
This indicates that two pairs of solutions bifurcate at the same value of A. This is partly a trivial 
observation, however, for the symmetry in the problem shows that solutions other than those for 
which u (|) = occur in symmetric pairs. The essential nature of the bifurcation can be seen 
by looking only at the left branch of the bifurcation curve shown in Figure 9. 

For comparison we also show a graph, Figure 10, of G for the forcing function g (t) = cost — 
sin2t. This is computed with A a bit above the bifurcation point. In this problem symmetry 
is lost, and solutions with u' = at and tc are not necessarily periodic. However the problem 
of finding solutions satisfying these two boundary conditions is still of interest, as an example of 
steady-states for the corresponding pde boundary value problem. The figure suggests that for 
some g the bifurcation may be of pitchfork type. We used the program xpp of G. B. Ermentrout 
ERMfl. 
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Figure 9 




Figure 10 

To study the bifurcation analytically we again use v — We are considering only the 
specific equation ([!]) . In order to prove that the bifurcation is not pitchfork, we study the 
linearization around the antisymmetric solution u p . We determine the slope of the function G 
at the antisymmetric solution. Since G (a) = u' a (ir), G' (a) = v' a (ir). But in the proof of 
Theorem we showed that when u = u p , v > on [0, 7r] and v' (tt) > 0. This proves that 
the bifurcation is not of pitchfork type, because the new solutions must appear at a positive 
distance away from u p . 

However, this analysis does not eliminate the possibility of several bifurcation points, in which 
solutions appear and then disappear. To get a complete picture more work is required. We have 
seen that for A > Ao there are at least five solutions. We also showed that At, the bifurcation 
point, is less than A . But in the following result, we hold A fixed and less than Ao, and let e tend 
to zero. 

Theorem 37 Suppose that A < Ao- Then for sufficiently small the only solution with 

u' (0) = 0, u' (tt) = 0. 

Corollary 38 lim e ^ \ b = A . 

Proof. The corollary follows immediately from the Theorem, the proof of which is more 
easily understood by reference to figure 1. We use the original scaling ([!]) . For any 5 > there 
are fi > and v > such that / (t, u) < —v on the set 

n Stfl = {(t,u) | (1 -//)?r < t < (1 + fi)n, u<U (t) - 5} . 

Hence, for sufficiently small e, no periodic solution can intersect the region Qs,tJ, at a point 
where (l — n < t < ( 1 + 7r, since then u would be forced below —b before \t — ir\ = fi. 
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So any periodic solution must have lie above U (t) — 5 on (l — j) tt < t < (l + |) it. The 



argument used to prove Theorem [M] can then be used to show that for sufficiently small e, u p is 
the only periodic solution u a with a maximum at tt. This completes the proof of the theorem 
and Corollary. ■ 

We now show that for small e there is no "reverse bifurcation" as A increases from A&. In 
other words, as we increase A from A&, there is no return to the case where u p is the unique 
periodic solution. 

Theorem 39 For sufficiently small e > 0, if X > A& then there are at least five periodic solutions, 
Mi, us with a, = Ui (0) < 0. Further, Ui < u.- l+ i for i = 1, • • • , 4. Here U3 is the antisymmetric 
solution, also denoted by u p . The solutions u\ and M5 are symmetric reflections, with w 5 (t) = 
—Mi (t + it) . Similarly, w 4 (t) = —u 2 (t + tt). 

Proof. We have seen that G (a) < for large negative a and also G < just below a p . The 
existence of at least five solutions for A > A& follows from the remarks above by showing that if 

«i = inf {a I u' a (it) = 0} , 

and if G (ax) = G' (at\) = 0, then 9G q^ > 0. Let u = u ai , h(t) = ^^-\ a=ai , and as above, 
v = tt. We see that 

e 2 h" = (3m 2 -X)h-u (36) 
h (0) = ti (0) = 



while v satisfies Q3"3]). Therefore we obtain, for any t G [0, 



t 

hv' — vh'\t = — I v (s) u (s) ds. (37) 



Lemma 40 For sufficiently small e, u ai < on [0, it] . 

Proof. This follows by a slight modification of the proof of Theorem |37|. Let 

(l StlM = {(t, u) I (1 - n) it < t < (1 + fjt) it, < u < U (t) - 8} . 

Suppose A = A . Then for small enough 5 and u, there is an e such that if < e < e , then any 
solution which intersects must decrease monotonically, in at least one direction, to below —6, 
within a time which is bounded over fls tf i- By continuity, this is also true for A sufficiently close 



to Aq. As in the proof of Theorem 0, no solution other than u p can have a positive maximum 



at 7T and not intersect Clg^, proving the result. ■ 

Also, if u = u ai , then v' (it) = G' (cui) = 0. Therefore v (it) 7^ 0. Suppose that v (it) < 0. 
Then for a slightly lower than a\, u a (tt) > u ai (tt) . Let 

P = inf {a I u a (t) > u ai (t) for some t G (0, 7r]} . 
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Then (3 is well-defined and — oo < (3 < a.\. Suppose that up (to) = u ai (to) for some to £ (0, it) . 
Since different solutions cannot be tangent, we must have up (t) > u ai (t) for some t G (0, tt) , 
but this contradicts the definition of f3. Therefore Up (tt) = u ai (tt) and u'^ (tt) > u ai (tt) = 0. 
But in this case, we can lower a further, until we find a 7 < (3 with u'^ (tt) = 0. This contradicts 
the definition of a±. A similar argument shows that v > on [0, tt]. 

Therefore, when u = u ai , v (tt) > 0. Then ([37]) and Lemma ^ show that du g ^ = h! (tt) > 0. 
Thus, for A just above A& there are at least five solutions (using symmetry). Further, there 
cannot be a decrease to fewer than five as A increases further, for at any point where G = G' = 
we would again get h' (n) > 0. 



To complete the proof of Theorem [39] it is convenient to let ct 2 = sup {a < a p \u' a (tt) = 0} . 
(Numerically, it appears there is only one 27r-periodic solution between u\ and u p .) Then reflection 
and translation of ii Ql and u a2 by the transformation u (t) — ► — u (t + tt) give the additional two 
asymmetric solutions. 

Our construction implies that Ui (0) < Ui+i (0) . From Proposition ^| and the way we define 
«2 we have u' a (tt) < for a 2 < a < a p , and in this range and close enough to a Pi u a < u v 
(This is because when u — u p> v > on [0, tt] .) Suppose, however, that for some a G [a2,cn p ), 
u a (t) > u p (t) for some t G [0, tt] . Let a = sup {a G [a 2 , a p )\u a (t) = u p (t) for some t G [0, tt].}. 
By the same argument as above we show that u& (tt) = u p (tt) and u' & (tt) > 0. Hence there is 
an a G (a, a p ) with u' a (tt) = 0. But this contradicts the definition of 

This proves that u 2 < u p . The proof that u\ < u 2 is similar, and the construction of U4 and 
«5 by reflection and translation implies the remaining order relations, namely, u p = n 3 < m 4 < it 5 . 
This completes the proof of Theorem |3^. ■ 



6 Conclusion 



6.1 How special is the cosine? 



General non-symmetric forcing functions are a topic for further study, but a few things are 



easy to see. First, the procedure in Theorems [T^ and [15] for obtaining some sort of chaotic 
behavior will carry over to a large variety of forcing functions. Also, on a finite interval the 
technique will give many steady states for the problem (0) — (^) . The only requirement is 
the existence of some set of functions Wk whose graphs form "fingers" pointing up and down 
alternately in a way similar to that in figure 2. We will not try to formulate a precise result 
here. 

We did use symmetry essentially to obtain the existence of the solution u p . It is here that we 
have found significant differences between the problem of finding steady-states for (|^) - (|) 011 a 
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finite interval and the problem of finding bounded solutions, and chaos, on an infinite interval. It 
is easy to find a third solution to the boundary value problem on a finite interval using shooting. 
From there one can go on to find many other solutions, both stable and unstable. But dealing 
with chaos on an infinite interval seems different, and it is for that reason that we concentrated 
on the particular equation (H). Having done so, it is natural to make use of symmetry to obtain 
simpler proofs in some cases where a shooting method may apply even without symmetry. 



We note, however, that the proof of Theorem 12, which does not use symmetry, includes 



the existence of a solution which does not intersect any of the w^. This solution is the stable 



solution of | AMPP |, and could play the role of u p in a study which uses shooting but does not 



use symmetry. 

The proof of stability of the solution u p (Theorem ^4[) does not use symmetry, but does use 
the fact that the derivative of cosine is negative in (0,7r). We do not believe that this property 
should be crucial to obtaining an elementary proof of linearized stability, and we expect that our 
proof will extend to more general forcing functions, but this is the subject of future investigation. 
The assumptions in | |AMPP|| imply that the forcing function g is monotone in a neighborhood of 



any of its zeros, and this should allow our proof to go through. However the result will probably 
be weaker in that Sg will be restricted more. 



We make use of Lemma about the monotonicity of the maxima and minima in (0, 7r), a 
number of times, and this also relies on the monotonicity of cosine in this interval. Our future 
work will include a study of how results will change without this condition, but we expect that 
extensions will be possible. 

The results about bifurcation as A increases, however, may change radically if g has less 
symmetry. For example, if g has period 2n, and is symmetric around ir, but not anti-symmetric 
around |, then the solutions will not appear in symmetric pairs. The technique we have used 
may still be able to prove that saddle-node bifurcations occur, at least in a nearly-symmetric 
situation, but there may be two of them, producing first two new solutions and then two more. 
This is easily seen in numerical simulations. As the deviation from full symmetry increases, 
it appears from some brief numerical experiments that other possibilities exist, and we hope to 
explore this further. 



6.2 Summary of main points. 

In Section 1 we introduce the problem and relate it to some previous work. 

In Section 2 we give preliminary results valid for all e > 0. Some of these are for a more 
general forcing function g (t) . The main points are that for A < there is a unique bounded 
solution, which is periodic, while for A > Ao there are at least three solutions. Therefore a 
bifurcation takes place, and some computations are given in section 4.5 indicating that this can 
be of different types depending on g (t) . 
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In section 3 we give results valid for a specific range of e and A. The main hypothesis is 



Condition [TTJ, the existence of "spikes", ( the w^), which are solutions tending to ±00 in both 
directions. Solutions are characterized by which of the Wk they intersect. Section 3.1 gives 
the proof of the 1:1 correspondence with certain sequences if Condition [ll| holds. It would be 
possible to rephrase this result to give a natural correspondence with sequences of three symbols, 
corresponding to the three solutions which we later labeled U\, u 3 (= u p ) , and w 5 . In 3.2 we give 
a brief discussion of "kneading theory" in our context. In Section 3.3 the symbolic dynamics is 



extended to sequences of five symbols. In section 3.4 Condition [T_TJ is verified, first for "sufficiently 
small" e, where no analysis is required, and then for larger e. 

The results in Section 3 do not include uniqueness or stability, and there is only a limited 
sensitivity to initial conditions demonstrated. They do not depend at all on symmetry, and 
indeed, the techniques will yield a weak form of chaos for a wide variety of forcing functions, 
including non-periodic forcing. 

Section 4 contains a variety of results, all proved for sufficiently small e with no estimate on 
the range of e for which they hold. In 4.1 some results are given about asymptotic behavior of 
solutions as e —>■ 0. In 4.2 further periodic solutions are found, including solutions with multiple 
internal layers and with a "down-jump" near |, in contrast to the solution u p which jumps 
upward at |. In 4.3 we give a proof of the stability of u p using classical ode methods, and 



extend the result on uniqueness of ||AMPP|| a bit by obtaining a larger region in which it is 



unique. The uniqueness proof is also more direct than that in 

In 4.4 we discuss sensitivity with respect to initial conditions. In 4.5 we consider the bifur- 
cation problem in A for 27r-periodic solutions of ([[]). 

Finally, in Section 5.1, we discuss the role of the specific cosine forcing function in our results, 
and conjecture that its symmetry and monotonicity properties may not be essential except in 
the bifurcation analysis of section 4.6 . 

We wish to express our appreciation to Professor H. Matano for sharing his thoughts on this 



problem. His remarks are cited in more detail at the end of section |5.2| , together with a citation 
of recent work of his coworker, Dr. K. Nakashima. 
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